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A  TRANSPORT  AND  DIFFUSION  MODEL  FOR  SMOKE  MUNITIONS 


1.  INTRODUCTION 

The  AMSAA  smoke  model ,  referred  to  as  the  Smoke  Effec¬ 
tiveness  Manual  Model  (SEMM),  as  described  here  was  derived  from 
the  original  JTCG  Smoke  Obscuration  Model  No.  1  (SOM-1)  whose  principal 
author  was  M.  C.  Johnson.  In  order  to  make  the  model  more  manageable, 

Mr.  G.  Hanna  of  AMSAA  and  the  Oklahoma  State  University  Engineering 
Field  Office  at  Eglin  Field,  Florida,  revised  the  SOM-1  by  removing  the 
more  doubtful  transmission  and  brightness  calculations  leaving  this 
information  to  be  supplied  by  calibrated  test  results.  Also,  additional 
features  were  included  such  as  the  effects  of  weapon  delivery  errors  and 
a  comprehensive  analysis  of  the  density  distribution  of  the  smoke  cloud. 

The  AMSAA  smoke  model  forms  the  basis  for  constructing  the 
obscuration  tables  given  in  the  Smoke  Effectiveness  Manual  published 
by  the  Joint  Technical  Coordinating  Group  for  Munitions  Effectiveness. 

It  considers  the  delivery  of  white  phosphorus  (WP)  or  hexachl orethane 
(HC)  muntions  by  indirect  fire  weapons  to  selected  aimpoints  located  at 
a  given  range  from  the  delivery  weapons.  Single  or  multiple  volleys 
may  be  fired.  After  impact,  the  computer  program  determines  the  amount 
of  obscurant  at  various  time  intervals  along  numerous  lines  of  sight. 

The  obscuring  screen  is  transported  and  diffused  as  a  function  of  local 
meteorological  conditions  during  which  time  a  criterion  is  applied  to 
determine  if  target  detection  can  be  achieved  with  the  particular  visual 
aid  employed.  The  model  is  adaptable  to  a  number  of  spectral  ranges  and 
visual  aids  depending  on  the  input  data  used.  The  following  distinctive 
features  are  noteworthy: 

•  The  smoke  model  is  a  transport  and  diffusion  model  and  re¬ 
quires  transmission  data. 

•  The  model  assumes  an  uncorrelated  guassian  tri variate 
distribution  for  each  obscuring  burst. 

•  The  model  produces  "holes"  or  discontinuities  in  the  smoke 
screen  due  to  the  aiming  and  precision  errors  of  indirect 
firing  weapons. 

•  The  model  is  used  for  detection  but  with  the  proper  data 
can  be  used  for  recognition  and  identification. 

According  to  Johnson,  Reference  2,  the  gaussian  trivariate 
distribution  of  aerosols,  chosen  for  this  model,  evolved  from  the  works 
of  Sutton,  Cal der  and  Milly.  Milly  compared  the  model  against  numerous 
test  data  on  munitions  of  various  types  in  order  to  establish  the  values 
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of  the  model's  constants.  Johnson  states  that  the  model  is  usually 
selected  by  analysts  for  short  to  medium  cloud  travels  because  of  the 
simplicity  of  parameters  required  for  input  information. 

A  recent  adjustment  has  been  made  to  the  WP  portion  of  the 
model  and  is  currently  being  evaluated.  This  change  will  improve  the 
correlation  of  the  model  output  with  field  data  from  smoke  tests. 
Documentation  will  be  published  at  a  later  date. 


2.  MODEL  DESCRIPTION 
2.1  General . 

The  AMSAA  model  can  be  used  for  both  WP  and  HC  munitions. 

Some  model  changes  are  required  when  employing  either  agent  but  a  basic 
similarity  exists  between  both  applications. 

When  bursted,  a  WP  munition  develops  an  initial  size  which 
proceeds  to  diffuse  and  transport  according  to  meteorological  conditions. 
Similarly,  the  smoke  produced  by  HC  munition  is  treated  as  a  large  num¬ 
ber  of  small  individual  continuous  bursts  or  puffs  which  form  an  elon¬ 
gated  screen  with  a  fixed  origin.  In  both  cases,  the  bursts  are  assumed 
to  have  a  tri variate  gaussian  distribution  with  standard  deviations  a^, 
Oy,  Ojr  at  any  time  t  given  by  Sutton  as: 
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=  0.1522  ( 


Ut  +  A^.9294 

1.0  ^ 
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=  3.41  ( 
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(2.1) 


The  quantity  U  represents  the  mean  wind  velocity  near  the  ground  and 
t  measures  the  age  of  each  burst.  The  exponents  a  and  3  are  func¬ 
tions  of  the  temperature  gradient  in  the  vicinity  of  the  ground.  The 
quantities  A,  B,  and  C  account  for  the  size  of  the  initial  burst  of  the 
munition.  These  values  are  zero  for  HC  rounds  but  have  the  following 
values  for  WP  rounds: 
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where  Oxs*  °ys»  °zs  s'"®  the  source  values  of  the  standard  devia¬ 
tions  and  depend  on  the  WP  fill  weight  of  the  munition.  The  constants 
in  Equations  (2.1)  and  (2.2)  were  established  by  Milly,  Reference  1,  in 
accordance  with  test  data  but  have  been  modified  somewhat  by  past  users. 
The  growth  of  these  sigmas  with  time  establishes  the  diffusion  qualities 
of  the  bursts  and  consequently  that  of  the  smoke  screen. 

For  calculating  convenience,  the  model  employs  three  coordi¬ 
nate  systems;  these  are  (a)  an  earth  fixed  system,  (b)  a  weapon  axis 
system  with  the  y  axis  extended  in  the  direction  of  the  line  of  fire, 
and  (c)  an  axis  system  so  oriented  that  its  x  axis  is  always  along  the 
wind  direction.  The  latter  facilitates  the  integration  of  smoke  density 
along  the  selected  lines  of  sight. 

The  burst  points  of  the  munitions  differ  from  the  intended  aim 
points  by  the  aiming  and  precision  errors  in  firing.  Figure  1  shows  the 
intended  deployment  of  munitions  in'  a  lazy  W  pattern  along  with  the 
actual  bursts  as  delivered  by  a  battery  of  six  weapons,  firing  at  2/3 
their  maximum  range,  on  one  occasion.  An  estimate  of  the  range  and 
deflection  deviations  of  the  actual  centroid  of  burst  points  from 
the  intended  centroid  of  bursts  due  to  the  aiming  error  is  estimated  by 
sampling  the  random  variable  x  from  a  uniform  source  and  proceeding  as 
shown  in  Figure  2  to  arrive  at  the  corresponding  miss  distances  R  and  D. 
These  distances  are  laid  out  along  the  weapon  system  axis. 

A  similar  procedure  is  followed  to  find  the  contribution  of 
the  precision  error  at  each  burst  point.  Components  of  the  precision 
and  aiming  error  are  summed  to  arrive  at  the  total  error  between 
the  intended  and  actual  burst  points.  Range  and  deflection  standard 
deviations  for  both  errors  must  be  known. 

The  inclusion  of  munition  delivery  errors  is  essential  to  the 
development  of  "holes"  or  periods  of  relatively  short  visibility  through 
the  smoke  screen.  Precision  errors  are  mainly  instrumental  in  creating 
"holes"  in  predomi nently  cross  wind  cases  whereas  aiming  errors  can  pro¬ 
duce  significant  smoke  screen  displacements  in  obscuring  an  enemy  front 
in  prevailing  head  or  tail  wind  situations. 

During  the  development  of  the  smoke  screen  from  the  individual 
bursts,  the  density  of  the  screen  at  any  point  is  the  sum  of  the  densi¬ 
ties  of  the  individual  burst  taken. at  that  point.  With  a  mean  ground 
wind  velocity  U  in  a  direction  x,  the  cloud  density  at  point  x,  y,  z 
in  the  wind  axes  system  is  given  by: 
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FIGURE  1  -  MUNITION  DISPERSION  DUE  TO  PRECISION  AND  AIMING  ERRORS 
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=  range  std  deviation 


B  =  Zttx 
R  =  A  sin  Bof^ 

D  =  A  cos  Bcjp 

FIGURE  2  -  DEVELOPMENT  OF  PRECISION  AND  AIMING  ERRORS 
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where  X  is  the  munition  efficiency  with  which  the  smoke  producing 
material  is  used. 


Q  is  the  weight  of  the  smoke  producing  material  prior  to 
dissemination,  excluding  weight  lost  in  the  plume. 


n  is  the  yield  factor  associated  with  the  physio-chemical 
reaction  process  which  converts  smoke  producing  material 
into  smoke.  For  hydroscopic  agents  this  quantity  is 
mainly  a  function  of  relative  humidity. 

z(t)  is  the  function  which  describes  the  vertical  motion 
of  the  smoke  due  to  updrafts  from  heat  released  dur- 
production. 

N  is  the  number  of  bursts  of  WP  or  puffs  of  HC  forming 
the  smoke  screen. 


The  integration  of  the  density  Equation  (2.3)  is  carried  out 
along  various  lines  of  sight  L-j  originating  at  the  observer's  position 
and  passing  throught  the  screen.  The  mass  per  unit  area  of  the  ob¬ 
scurant  along  line  L^-  is  therefore: 


MASS  (L)  = 


(DEN)  dL 


(2.4) 


The  line  of  sight  1]  is  an  obscured  path  if  and  when: 


MASS  (Li )  >  CLTHRS 


(2.5) 


where  CLTHRS  is  the  threshold  value  for  the  combined  smoke  agent  and 
optical  aid  employed.  That  is,  it  is  the  mass  per  unit  area  through 
which  the  probability  of  target  detection  is  derived  due  to  the  at¬ 
tenuation  of  the  obscurant. 
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The  obscuration  characteristics  of  the  smoke  cloud  is  deter¬ 
mined  by  exploring  various  lines  of  sight  (LOS)  through  the  developing 
screen.  Figure  3  shows  a  screen  formed  by  the  diffusion  of  three  bursts. 
An  enclosing  rectangle  is  given  by  the  four  sigma  values  of  those  bursts 
which  define  the  rectangle's  maximum  dimensions.  The  size  of  the 
rectangle  gives  the  limiting  locations  of  the  exploratory  LOS's  and 
assumes  that  these  LOS's  form  unobscured  paths.  Figure  4  expands  this 
idea  indicating  the  individual  LOS's  along  with  their  respective  separa¬ 
ting  lengths,  Ax.  The  screen  length  is  defined  by  the  first  and  last 
LOS  which  satisfies  the  condition  given  in  Equation  (2.5).  Those  interior 
Ax's  which  correspond  to  those  LOSs  which  fail  to  meet  this  criterion 
are  summed  to  give  the  effective  size  of  the  "hole"  in  the  screen.  The 
effective  cloud  length  is  given  by  the  total  length  less  the  "hole" 
size.  Replications  of  the  above  procedure,  repeatedly  sampling  the 
precision  and  aiming  errors,  provide  the  average  values  of  the  total  and 
effective  cloud  lengths. 

The  model  makes  provisions  to  explore  the  characteristics  of 
the  smoke  cloud  in  greater  depth.  Thus  the  densities  at  numerous  points 
within  the  smoke  cloud,  at  a  given  height  above  ground  level,  are  pre¬ 
sented  as  an  array  of  data  thereby  making  it  possible  to  follow  the  time 
history  of  high  and  low  density  areas.  Another  feature  is  the  calcula¬ 
tion  of  the  smoke  concentration  along  numerous  lines  of  sight  originating 
at  points  along  the  near,  or  observer  edge,  of  the  cloud  and  terminating 
at  points  along  the  far  edge.  In  this  way  multi  observer-target  LOS's 
are  explored  for  each  cloud  time  history. 


3.  MODEL  INPUT  DATA 

The  model  inputs  needed  to  produce  obscuration  results  may  be 
divided  into  two  parts,  selected  inputs  and  fixed  inputs.  The  selected 
inputs  are  those  which  the  user  can  choose  in  battlefield  planning  in 
order  to  employ  the  JTCG/ME  Smoke  Manual.  They  are  also  the  most  signifi¬ 
cant  parameters  determining  the  history  of  the  smoke  cloud.  The  fixed 
inputs  are  those,  which  have  a  lesser  influence  on  the  dissemination  of 
the  smoke  cloud  and  which  are  not  easily  discernible  or  calculable  on 
the  battle  field.  They  are  inputs  that  are  either  held  constant  in 
developing  the  smoke  tables  or  values  that  are  dependent  on  those  of  the 
selected  inputs. 

3.1  Selected  Inputs. 

3.1.1  Meteorological  Conditions.  Select  meteorological 
conditions  have  been  reduced  to  five  combinations  of  atmospheric  stability 
and  wind  velocities  in  accordance  with  the  Pasquill  grouping  of  these 
parameters.  These  are  in  Table  1. 
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FIGURE  4  -  LINE  UP  SIGHT  EXPLORATORY  REGION 


TABLE  1.  METEOROLOGICAL  CONDITIONS 


Stability 

Wind  Speed, 
knots 

Approx. 

Pasquill 

Category 

Ideal  (Inversion) 

5 

E  to  F 

Favorable  (Neutral) 

5 

D  to  C 

Favorable  (Neutral) 

10 

D  to  C 

Favorable  (Neutral) 

15 

D  to  C 

Marginal  (Lapse) 

5 

A  to  B 

3.1.2  Wind  Direction.  The  wind  directions  are  chosen  with 
reference  to  the  observer-to-target  LOS.  These  are  crosswind,  head  or 
tail  wind,  and  quartering  wind. 

3.1.3  Munition  Type.  The  model  will  produce  obscuration  re¬ 
sults  for  currently  available  WP  and  HC  munitions.  Those  to  be  included 
in  the  Smoke  Manual  are  shown  in  Table  2. 

TABLE  2.  SELECTED  SMOKE  MUNITIONS 


105mm,  M84B1,  HC 
105mm,  M60A2,  WP 
155mm,  M110E2,  WP 
155mm,  M116B1,  HC 
4.2",  M328A1,  WP 
81mm,  M375A2,  WP 
5",  5.54  CAL.,  WP 

3.1 .4  Vol ley  Size.  One  or  more  multi- round  volleys  may  be 
fired  at  selected  times  between  volleys.  To  avoid  lengthy  tables  in  the 
Smoke  Manual,  only  the  effects  of  single  round  or  multi-round  volleys 
will  be  listed  and  the  effects  of  additional  volleys  approximated  by 
simple  manual  calculations. 
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3.1.5  Spectral  Bands.  The  model  is  adaptable  to  any  spectral 
band  for  which  the  sensor's  capability  against  a  particular  smoke  agent 
is  available.  This  capability  is  represented  by  the  concentration- 
length  of  smoke  which  the  sensor  can  defeat.  For  the  Smoke  Manual, 
results  will  be  listed  for  the  visible  spectrum,  for  an  anti-tank  guided 
missile  and  for  the  currently  available  night  sight  operating  against 
the  smoke  munitions. 

3.2  Fixed  Inputs. 

As  stated  in  paragraph  3,  the  fixed  inputs  are  those  which  are 
either  held  constant  in  preparing  the  smoke  tables  or  whose  values  are 
average  values  corresponding  to  the  range  of  selected  inputs. 

3.2.1  Yield  Factor.  Most  military  screening  smokes  consist 
chiefly  of  hydroscopic  products  of  chemical  reactions  which  form  dilute 
solutions  in  the  presence  of  atmospheric  water  vapor.  The  resulting 
smoke  clouds  are  actually  suspensions  of  small  liquid  droplets.  Both  WP 
and  HC  munitions  are  in  this  class  and  the  mass  of  smoke  they  produce  is 
greater  than  the  mass  of  the  original  dry  agent.  Fog  oil  smoke,  on  the 
other  hand,  is  an  atomization  of  a  liquid  and  the  smoke  mass  is  equal  to 
the  mass  of  oil  consumed.  The  yield  factor  is  the  ratio  of  the  wet 
smoke  mass  produced  to  the  mass  of  dry  agent  consumed.  Yield  factors 
for  various  smoke  agents,  taken  from  Reference  2,  are  shown  in  Figure  5 
plotted  as  a  function  of  relative  humidity.  A  constant  relative 
humidity  of  50  percent  was  used  in  developing  the  smoke  tables. 

Table  3  gives  values  of  yeild  factor  used  in  developing  the  smoke  tables. 


TABLE  3.  YIELD  FACTORS  FOR  50  PERCENT  RELATIVE  HUMIDITY 


Type  Munition 

WP 

HE 


Yield  Factor 

5.2 

1.85 


3.2.2  Diffusion  Parameters  a  and  g.  According  to  Milly, 
Reference  1,  the  diffusion  of  the  smoke  bursts  (or  puffs  in  the  case  of 
HC)  are  controlled  by  the  meteorologically  dependent  quantities  a  and 
3  which  are  determined  by  the  local  temperature  gradient.  A  reproduc¬ 
tion  of  Milly's  curve  is  shown  in  Figure  6.  a  and  3  develop  the 
standard  deviations  for  the  trivariate  smoke  density  distributions  in 
directions  perpendicular  to  the  wind  direction.  Along  the  wind  direc¬ 
tions,  the  diffusion  parameter  6  has  the  constant  value  0.9294.  The 
specific  temperature  gradients  and  corresponding  a  and  3  values  are 
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FIGURE  5  -  YIELD  FACTOR  AS  A  FUNCTION  OF  RELATIVE  HUMIDITY  FOR 
VARIOUS  SMOKE  PRODUCING  SUBSTANCE'S 


a  =  10 
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FIGURE  6  -  (a  =  3)  vs  TEMPERATURE  DIFFERENCE 


given  in  Table  4  for  the  three  atmospheric  stability  conditions  covered 
in  the  smoke  tables.  The  gradient  is  represented  by  the  average  tempera¬ 
ture  difference  between  a  height  of  4  meters  and  0.5  meters  along  the 
local  terrain. 


TABLE  4. 

a  AND  B  VALUES  FOR  THREE 

STABILITY  CONDITIONS 

AT.  op 

a  =  B 

Stability  Conditions 

-  1.5 

1.5 

Lapse 

0 

0.88 

Neutral 

+  1.5 

0.70 

Inversion 

The  stability  conditions  are  selected  inputs. 

3.2.3  Munition  Fill  Weight.  The  AMSAA  Model  uses  the  munition 
fill  weight  modified  by  an  efficiency  factor  to  account  for  two  effects; 
(1)  the  amount  of  fill  actually  converted  to  dry  smoke  and  (2)  the  amount 
of  smoke  which  is  lost  in  the  rising  plume.  The  former  is  determined 
from  test  results  while  the  latter  is  dependent  on  atmospheric  stability 
and  is  employed  instead  of  using  a  smoke  rise  function.  The  product  of 
fill  weight  and  efficiency  is  called  the  effective  fill  weight  for  the 
munition  and  its  values  are  given  in  Table  5  for  the  munitions  listed 
in  Section  3.1 .3. 

TABLE  5.  EFFICIENCIES  AND  FILL  WEIGHTS  FOR  VARIOUS  MUNITIONS 


Efficiency  Effective  Fill  Wt.  (lb) 

Fill  Wt. 


Munition 

(lb) 

Lapse 

Neutral 

Inversion 

Lapse 

Neutral 

Inversion 

105  WP 

3.85 

0.10 

0.10 

0.30 

0.38 

0.38 

1.15 

105  HC 

7.50 

0.75 

0.75 

0.75 

5.60 

5.60 

5.60 

155  WP 

15.6 

0.10 

0.10 

0.30 

1.56 

1.56 

4.68 

155  HC 

25.8 

0.75 

0.75 

0.75 

19.4 

19.4 

19.4 

4.2  in. 

7.5 

0.40 

0.40 

0.60 

3.0 

3.0 

4.5 

WP 

81mm  WP 

1.75 

0.10 

0.10 

0.30 

0.175 

0.175 

0.53 

5  in.  WP 

8.34 

_ 
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3.2.4  Source  Sigmas.  °yci  ^yg.  °7c.  While  HC  smoke 
originates  from  a  small  area  which  may  be  considered  a  point  source,  WP 
smoke  originates  as  the  result  of  a  bursting  charge  which  ignites  and 
separates  the  filler  material.  The  distribution  of  this  material  is 
assumed  to  have  the  same  gaussian  characteristics  as  the  subsequently 
diffused  cloud.  The  source  sigmas,  defined  as  ®xs»  ®ys»  ®zs 
Equation  (2.2),  characterize  the  initial  WP  burst  and  have  values 
strongly  dependent  on  the  munition's  effective  fill  weight.  Their 
growth  with  time  defines  the  diffusion  process  of  the  smoke  screen. 

Under  calm  wind  conditions  ®xs  equal  to  ‘^yc  so  that  the  cloud 
obtains  a  circular  symmetry  in  planes  parallel  to  the  ground.  In  strong 
winds,  the  down  wind  sigma,  ^xs  skewed  in  the  wind  direction. 

For  HC  munitions  the  source  sigmas  are  given  a  value  of  zero.  Figure  7, 
taken  from  Reference  2,  shows  the  calm  wind  sigmas  plotted  against  the 
munition's  effective  fill  weight  as  determined  by  limited  test  results. 
For  the  WP  munitions  listed  in  Section  3.2.3,  the  source  sigmas  were 
used  in  preparing  the  Smoke  tables.' 


TABLE  6. 

SOURCE 

SIGMAS  FOR 

VARIOUS  WP 

MUNITIONS 

a,,  = 

®vc.  meters 

®7c.  meters 

Munition 

Lapse 

Neutral 

Inver. 

Lapse 

Neutral 

Inver. 

105  WP 

2.79 

2.79 

3.79 

0.93 

0.93 

1.26 

155  WP 

4.13 

4.13 

5.62 

1.37 

1.37 

1.87 

4.2  in.  WP 

4.96 

4.96 

5.56 

1.65 

1.65 

1.87 

81mm  WP 

2.24 

2.24 

3.05 

0.74 

0.74 

1.01 

5  in.  WP 

_ 

3.2.5  Light  Attenuation  and  Extinction  Coefficient.  The  frac¬ 
tional  loss  of  light  intensity  and  the  extinction  coefficient  are  assumed 
to  be  adequately  related  by  the  Beer-Lambert  Law.  When  the  threshold 
level  of  attenuation  and  the  extinction  coefficient  for  a  particular 
smoke  agent  are  used  in  this  relationship,  the  threshold  value  of  smoke 
concentration  times  length  is  obtained.  This  value,  noted  as  CLTHRS  in 
the  model,  is  the  model's  obscuration  criterion  and  is  determined  as  a 
function  of  the  selected  inputs.  In  fact,  the  ability  of  the  sensor  to 
"see"  through  the  cloud  of  a  particular  obscurant  depends  on  its  ability 
to  defeat  its  corresponding  concentration  -  length  threshold.  The  fixed 
inputs  of  attenuation  threshold  and  extinction  coefficients  used  in  the 
model  are  given  in  Table  7  for  the  visual  range,  for  an  anti-tank  guided 
missile  (ATGM)  and  for  an  operating  set  of  night  sights. 
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FIGURE  7  -  SOURCE  SIGMAS  VS  FILL  WEIGHT 


TABLE  7.  ATTENUATION  THRESHOLD  AND  EXTINCTION  COEFFICIENTS 


VISIBLE 

ATGM 

NIGHT  SIGHT 

ATTEN.  THRES 

• 

0.90 

0.99 

0.95 

WP 

0.00332 

0.00267 

0.00028 

HC 

0.0045 

0.0029 

0.000061 

3.2.6 

Munition  Delivery  Accuracy.  The  smoke  tables  are  con- 

structed  for  munition  delivery  accuracies  corresponding  to  two-thirds  of 
the  delivery  weapon's  maximum  range.  Representative  values  of  the  aiming 
and  precision  errors  for  the  cannon  and  mortar  weapons  covered  in  the 
manual  are  listed  in  Table  8.  The  aiming  errors  represent  adjusted  fire. 

TABLE  8. 

AIMING  AND 

PRECISION  ERRORS 

Range  (m) 

Precision,  Std.  Dev. 

Aiming,  Std.  Dev. 

Range 

Defl. 

Range  Defl. 

105  How. 

7000 

23 

6 

36  27 

155  How. 

12000 

44 

10 

46  15 

4.2  In.  Mort. 

4000 

43 

5 

49  27 

81mm  Mort. 

5000 

37 

6 

45  27 

5  In.  Gun 

16000 

3.2.7  Aimpoints.  The  intended  center  of  impact  for  the 
volley  located  from  the  center  of  the  firing  weapons  as  the  origin  of 
coordinates,  is  designated  as  YT  in  the  range  direction  and  XT  in  the 
lateral  direction.  The  values  for  YT  correspond  to  the  weapon  ranges 
quoted  in  Section  3.2.6  while  XT  is  taken  as  zero.  The  location  of  the 
intended  impact  points  about  the  intended  center  of  impact  taken  as  the 
origin,  are  noted  as  XIDEAL  or  X(I)  in  the  lateral  direction  and  YIDEAL 
or  Y(I)  in  the  range  direction.  These  points  also  establish  the  desired 
impact  pattern.  Table  9  gives  the  intended  impact  points  for  some 
of  the  volley  sizes  used  in  constructing  the  tables  in  the  Manual. 
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TABLE  9.  INTENDED  IMPACT  POINTS  FOR  VARIOUS  MUNITIONS  AND  VOLLEY  SIZES 


Firing 

Voll ey 

XIDEAL 

Y I DEAL 

Pattern 

Size 

(m) 

(m) 

105  How.  Lazy  W 

6 

-75 

-15 

-45 

15 

-15 

0 

15 

0 

45 

15 

75 

-15 

155  How.  Lazy  W 

6 

-125 

-25 

-  75 

25 

-  25 

0 

25 

0 

75 

25 

125 

-25 

4.2  In.  Line 

4 

-  60 

20 

Mort. 

-  20 

0 

20 

20 

60 

0 

81mm  Mort.  Line 

3 

-  30 

0 

0 

0 

30 

0 

5  In.  Gun  Converge 

2 

0 

0 

Sheaf 

0 

0 

3.2.8  Reliability.  The  reliability  of  all  munitions  is  assumed 

to  be  100  percent.  For  a  round 

rel iabil ity 

other  than 

100 

percent,  the 

model  would  produce  an  average 

cloud  length 

.  For  small 

volley  sizes,  say 

a  one  round  volley,  an  average 

value  would 

not  be  realized 

since  a  full 

cloud  either  would  or  would  not 

occur. 

3.2.9  Burn  Time.  The  burn  time  for  all  HC  rounds  is  assumed  to 
be  2  minute?!  This  is  in  accordance  with  field  test  data.  (Reference  4). 

4.  SAMPLE  RESULTS 


Sample  calculations  were  made  for  a  battery  of  155mm  weapons 
firing  in  a  cross  wind  of  5  knots  under  neutral  atmospheric  conditions. 
Table  10  shows  the  results  of  firing  a  single  volley  of  2,  4,  and  6  rounds 
of  WP  and  HC  ammunition.  The  battery  fire  has  been  adjusted  thereby 
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minimizing  the  aiming  errors.  The  cloud  characteristics  are  defined  as 
follows: 


(a)  Total  cloud  length  -  The  included  distance  between  the 
extreme  left  and  right  obscured  LOS  from  the  observer's  position. 

(b)  Effective  cloud  length  -  The  total  cloud  length  less  the 
sum  of  the  included  "hole"  lengths. 

(c)  Upwind  edge  of  cloud  -  The  distance  of  the  upwind  edge  of 
the  cloud  from  the  intended  volley  aimpoint  (centroid).  This  distance 
becomes  negative  when  the  upwind  edge  crosses  the  aimpoint. 

(d)  Downwind  edge  of  cloud  -  The  distance  of  the  downwind 
edge  of  the  cloud  from  the  intended  volley  aimpoint.  This  distance  is 
positive  when  the  downwind  edge  is  downwind  of  the  aimpoint. 

The  two- round  volleys  are  fired  at  aimpoints  3  and  4  of  the 
Lazy  W  shown  in  Figure  1  while  the  four-round  volleys  use  the  outer 
aimpoints  1,  2,  5  and  6.  These  aimpoints  were  chosen  to  illustrate  the 
results  obtained  using  close  and  separated  bursts  and  are  not  neces¬ 
sarily  those  to  be  used  in  actual  firings.  The  Uteral  spacing  between 
each  sequential  aimpoint  is  50  meters,  all  of  which  are  located  at  two- 
thirds  of  the  weapons'  maximum  range  from  the  firing  position. 

Referring  to  the  WP  two-round  volley  data  in  Table  10,  the 
closeness  of  the  aimpoints  results  in  a  cloud  with  an  initially  small 
hole  between  the  two  bursts.  This  is  indicated  by  the  difference  between 
the  total  and  effective  lengths.  The  total  obscuring  length  increases 
to  202  meters  in  3.5  minutes  after  which  the  cloud  thins  out  and  the 
obscuring  length  starts  to  shorten. 

The  150-meter  spacing  between  the  number  2  and  5  aimpoints  for 
the  four-round  volley  creates  a  larger  separation  for  the  corresponding 
bursts.  Results  in  Table  10  show  that  the  accumulated  hole  size  starts 
off  at  about  100  meters  in  size  at  30  seconds  after  impact  and  closes 
almost  completely  in  4  minutes.  Meanwhile,  the  effects  of  diffusion 
causes  the  ends  to  thin  out  so  that  at  3  minutes  after  round  impact  a 
shortening  of  the  effective  length  occurs. 

The  six-round  volley  increases  the  effects  of  the  four-round 
volley  by  the  addition  of  the  two  middle  weapons  of  the  battery.  Since 
the  outer  weapons  are  fired  in  both  cases,  a  similar  cloud  length  would 
be  expected  but  with  a  smaller  accumulated  hole  size  due  to  the  filling 
effect  of  weapons  3  and  4.  Table  10  shows  this  to  be  the  case,  the  hole 
size  being  only  a  few  meters  in  length  throughout  the  cloud's  time 
history.  The  end  point  displacements,  or  distances  covered  by  the  screens, 
are  expected  to  be  the  same  whether  four  or  six  guns  are  firing. 
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The  HC  burst  points  are  a  continuous  source  of  smoke  until 
burnout  occurs  at  the  end  of  2  minutes.  At  this  time  the  upwind  edge 
of  the  cloud  starts  to  move  downwind.  No  holes  are  indicated  after 
30  seconds  of  burn  time.  Since  the  four-  and  six-round  volleys  fire 
the  same  outer  weapons  of  the  battery,  little  or  no  advantage  is  gained 
by  firing  the  six  rounds  for  the  240  seconds  indicated  in  the  Table. 

NOTE:  Data  used  in  calculating  results  for  the  HC  agent  have  been 
updated  since  this  table  was  prepared  and  are  included  in 
Section  3  of  this  report. 


5.  COMPARISON  BETWEEN  MODEL  AND  TEST  RESULTS 

Between  September  and  November  1977,  Dugway  Proving  Ground 
personnel  tested  certain  US  inventory  smoke  munitions.  These  tests 
were  requested  by  AMSAA  and  funded  by  the  Office  of  the  Project 
Manager  for  Smoke  and  Obscurants.  Dugway 's  report  on  these  tests  is 
given  in  Reference  3.  These  tests  included  most  of  the  smoke  munitions 
listed  in  the  JMEM/ME  Smoke  Manual  in  addition  to  burn  time  data  for  HE 
munitions  taken  from  wind  tunnel  tests.  For  these  reasons  they  were  used 
as  a  basis  for  comparison  with  results  obtained  from  the  Smoke  Effective 
ness  Manual  Model.  Table  11  is  a  summary  of  the  munitions  and  the  test 
conditions  from  which  data  were  selected  for  the  comparative  analysis. 

5.1  Test  Data. 

In  order  to  cover  a  wide  range  of  volume  of  delivered  smoke, 
test  data  from  single-  and  multi-round  volleys  of  155mm  and  105mm  WP 
smoke  projectiles  were  selected  from  Reference  3.  The  choice  of  WP 
rather  than  HC  munitions  for  comparison  between  test  and  predicted 
results  represents  a  more  severe  test  of  simulation  due  to  the  defini¬ 
tive  nature  of  the  WP  burst  and  the  individual  burst  histories  of  trans¬ 
port  and  diffusion  which  the  model  must  consider.  Table  12  includes 
test  data  required  by  the  model  as  input  for  simulating  the  test  results. 
These  data  are  identified  in  the  table  which  also  lists  computer  instruc¬ 
tions  appropriate  for  handling  this  information.  Finally,  the  table 
includes  input  parameters  that  are  functions  of  the  input  conditions. 

The  latter  are  given  in  Figures  5,  6,  and  7  and  are  taken  from 
Reference  2. 

5.2  Model  Changes  for  Comparative  Analysis. 

A  minor  alteration  of  the  WP  program  of  the  smoke  model  was 
necessary  to  produce  results  for  comparison  with  test  data.  This 
change  is  shown  in  Appendix  D  Flow  charts  titled.  Subroutine  CALC, 
Modification  for  Comparative  Analysis.  Essentially,  this  change 
determines  the  value  of  concentration-length,  CL,  of  the  smoke  cloud 
along  a  fixed  observer-to-target  line-of-sight.  These  CL's  are  re¬ 
corded  versus  time  from  burst. 


26 


TABLE  11:  SUMMARY  OF  TESTS  FOR  MODEL  VALIDATION* 
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REF:  Inventory  Smoke  Munition  Test  (Phase  Ila)  Final  Test  Report,  Report  No.  DPG-FR-77-314, 
June  1978. 


TABLE  12.  MODEL  INPUT  PARAMETERS 
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TABLE  12.  (CONTINUED) 
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5.3  Comparisons  Between  Model  and  Test  Results. 

The  WP  portion  of  the  smoke  model  was  run  using  test  data 
obtained  from  the  Smoke  Inventory  Tests,  Reference  3.  In  order  to  get  a 
broad  view  of  the  model's  predictive  capabilities,  the  tests  selected 
for  comparative  analysis  consisted  of  those  employing  one  and  three 
rounds  of  105mm  WP  M60A2  munitions  and  three  and  six  rounds  of  155mm  WP 
M110E2  munitions.  Each  group  of  rounds  were  statically  fired  simul¬ 
taneously  in  a  linear  array  parallel  to  the  observer's  line  of  sight. 

In  this  way  the  smoke  cloud  passing  the  observer's  LOS  varied  considerably 
in  concentration-length  from  test  to  test  rather  than  in  screen  length. 

A  time  history  of  smoke  concentration  was  recorded  as  the  cloud  passed 
the  observer's  LOS  and  a  comparison  of  this  quantity  was  made  with  the 
corresponding  predicted  value  obtained  from  the  smoke  model.  Figure  8 
shows  the  predicted  and  test  results  when  the  smoke  from  one  round  of 
105mm  WP  is  released  and  crosses  the  observers  LOS.  The  concentration 
length  represents  the  mass  of  smoke  occupying  a  volume  one  square  meter 
in  cross  section  and  extending  along  the  observer's  LOS.  The  dashed  and 
solid  markers  indicate  the  level  of  concentrations  needed  to  defeat  the 
naked  eye  and  an  anti-tank  guided  missile  (ATOM)  sensor  respectively. 

The  corresponding  defeat  level  for  night  sight  is  above  the  concentration 
levels  shown.  Comparisons  between  model  and  test  results  at  the  threshold 
levels  indicated  is  quite  favorable.  A  characteristic  of  most  test  and 
model  comparisons  is  the  tail -off  of  smoke  that  persists  after  the  main 
body  of  smoke  has  passed  downwind.  This  effect  is  apparently  caused  by 
the  filler  material  containing  solid  chunks  which  do  not  form  the  initial 
flash  but  continue  to  burn  on  contact  with  the  ground.  This  effect  does 
not  significantly  degrade  the  fidelity  of  the  model  until  larger  volleys 
of  the  larger  155mm  munition  are  involved.  In  such  cases,  the  model 
gives  conservative  results.  Figure  9  shows  a  similar  presentation  for 
three  rounds  of  105mm  WP.  A  favorable  comparison  between  model  and 
test  results  is  also  apparent  at  the  visual  and  ATOM  thresholds.  Here 
again,  the  concentration  level  to  defeat  the  night  sight  is  not  reached. 
The  threshold  level  for  the  night  sight  is  finally  exceeded  with  three 
rounds  of  155mm  WP.  This  is  shown  in  Figure  10  along  with  the  visual 
and  ATOM  threshold  levels.  The  correlation  is  favorable,  with  the 
model  being  a  little  conservative  at  the  lower  concentrations  due  to  the 
tail-off.  Figure  11  shows  comparative  results  for  6  rounds  of  155mm 
WP.  In  this  particular  test  the  tail -off  persists  for  a  considerable 
length  of  time  after  the  maximum  smoke  concentration  level  has  passed. 

From  the  comparisons  made  to  date  between  model  and  test  re¬ 
sults,  it  appears  that  the  model  predicts  the  duration  of  obscuration 
for  the  three  sensors  quite  well.  Difficulties  arise  where  the  WP  model 
predicts  higher  peak  values  of  CL's  and  shorter  duration  times  at  lower 
cl's.  Both  effects  are  apparently  due  to  some  of  the  filler  material 
not  being  consumed  in  the  initial  burst  and  burning  on  the  ground  in  the 
manner  of  an  HC  munition.  A  model  change  has  been  prepared  to  correct 
this  deficiency.  This  change  will  be  documented  at  a  later  date. 
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FIGURE  8  -  PREDICTED  VS.  TEST  RESULTS,  ONE  ROUND  105mm,  WP 


TIME  SEC 

FIGURE  9  -  PREDICTED  VS  TEST  RESULTS, 
3  ROUNDS  105mm,  WP 
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TIME  SEC 

FIGURE  10  -  PREDICTED  VS  TEST  RESULTS 
3  rounds  155mm  WP 
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FIGURE  n  -  PREDICTED  VS  TEST  RESULTS, 
6  ROUNDS  155nin  WP 
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APPENDIX  A 


DEFINITION  OF  SYMBOLS 


The  input  quantities  required  to  operate  the  smoke  model  are 

listed  and  defined  below. 

ANG  The  angle  between  the  earth  axis  and  the  weapon  axes  system, 

degrees.  (Taken  as  zero  in  the  program  listing). 

CATTN  The  attenuation  coefficient  used  in  the  Beer-Lambert  Law. 
Function  of  the  smoke  agent  and  transmitted  wavelength. 

DL  Line  increment,  meters. 

DT  DT  is  the  time  between  computer  computations  or  steps, 

seconds. 

EFF  Defined  in  Section  3.3.2,  Munition  Fill  Weight. 

INSNAP  Time  between  printouts  (SNAP)  divided  by  time  between 
computations  (DT). 

NRPV  The  number  of  rounds  per  volley  to  be  fired. 

NS  The  total  number  of  computer  time  steps  in  the  program. 

NSAMM  The  number  of  "holes"  or  see  through  locations  within  the 

smoke  screen. 

NSAMP  Total  number  of  replications  for  averaging  results. 

NT,  MT  NT  is  the  total  number  of  times  the  cloud  history  is 

examined,  or  looked  at  each  time  noted  as  MT,  and  separated 
by  SNAP  seconds. 

NVOL  The  number  of  volleys  to  be  fired. 

QMUN  Defined  in  Section  3.2.3,  Munition  Fill  Weight. 

SNAP  Time  increment  between  printouts  or  "looks",  i.e. ,  between 
MT  and  MT+1 ,  seconds. 

THRES  The  percent  of  light  leaving  the  target  which  is  required 
to  accomplish  detection. 

U  The  wind  velocity  in  meters  per  second. 
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VT 

Time  between  volleys,  seconds. 

WVUl, 

WVU2 

The  direction  cosines  between  the  earth  axes  and  the  wind 
axes. 

Xob.yqb 

The  observer  location  in  the  earth  axes  system  (usually 
0,0)  meters. 

XIDEAL(J) 

YIDEAL(J) 

ZIDEAL(J) 

The  aimpoints  of  all  J  rounds  to  be  fired,  measured 
from  the  volley  aimpoint  or  desired  pattern  center. 

XT,  YT 

The  location  of  the  aimpoint  for  the  volley  in  the 
weapon  axes  system,  meters. 

YIELD 

Defined  in  Section  3.2.1. 

a,  0,  a 

The  exponents  defining  the  standard  deviations  of  the 
smoke  cloud.  Functions  of  temperature  gradient. 

m «  ‘^AD 

The  occasion-to-occasion  aiming  error  in  range  and 
deflection,  meters. 

OBR .  OBD 

The  round-to-round  precision  error  in  range  and  de¬ 
flection,  meters. 

°xs »  °sy 

The  source  sigmas,  meters.  Establish  the  initial  size 
of  the  smoke  burst;  a  function  of  fill  weight. 
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APPENDIX  B 


DESCRIPTIONS  OF  SUBROUTINES  FOR  WP  MUNITIONS 


B.l  GENERAL 

The  subroutines  for  the  WP  program  are  described  in  this 
Appendix  using  a  combination  of  algebraic  and  FORTRAN  notation.  The 
order  of  presentation  is  essentially  that  in  which  they  are  called  for 
in  the  MAIN  listing.  Flow  charts  for  the  WP  program  are  in  Appendix  D. 

B.2  WP  PROGRAM,  MAIN 

The  function  of  MAIN  is  to  call  sixteen  subroutines  which  make 
up  the  WP  smoke  model.  In  addition,  the  calculation  of  the  threshold 
value  of  smoke  CL  for  a  given  visual  aid  is  computed  here  in 
accordance  with  the  Beer-Lambert  Law.  Simply: 


dl/dCL  =  CATTN  x  I 

TRANS  =  1  -  ATTN  =  I/Iq  =  exp  (-CATTN  *  CL) 

ATTN  =  1  -  exp  (-CATTN  x  CL)  (6) 

where  I/Iq  is  the  ratio  of  the  transmitted  to  the  incident  radiation 
intensity. 

CATTN  is  the  extinction  coefficient  for  the  smoke  agent. 

TRANS,  ATTN  are  the  transmission  and  attenuation  of  radiation  through 
the  smoke,  respect ivity. 

B.3  SUBROUTINE  STCL 

Subroutine  STCL  calculates  the  virtual  offset  distances  A,  B, 
C  which  are  needed  to  define  the  initial  burst  sigmas  to  be  used  in 
Subroutine  FRMCLD.  It  also  establishes  the  weight  of  smoke  produced 
from  the  selected  smoke  agent  through  the  calculation  of  the  quantity 
FACTOR. 

B.4  SUBROUTINE  FRMCLD 

Subroutine  FRMCLD  computes  the  sigma  time  histories  for  all 
the  bursts  which  make  up  the  smoke  cloud.  The  increase  of  these  sigmas 
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with  time  is  the  diffusion  mechanism  which  creates  the  cloud  growth.  The 
quantity  FACT  is  the  smoke  density  written  in  the  form  necessary  for 
inclusion  in  the  gaussian  formulation. 

B.5  SUBROUTINE  CLEAR 

Subroutine  CLEAR  is  employed  to  set  to  zero  the  initial  values 
of  the  statistical  variables  which  are  to  be  evaluated  in  Subroutine 
VEVAL.  This  step  forms  part  of  the  Monte  Carlo  averaging  process  to  be 
preformed  in  Subroutine  CALPRT  which  follows  VEVAL.  The  quantities  to 
be  initially  set  to  zero  are: 


Cloud  length 

SSL 

Cloud  length,  squared 

SSLS 

Left  end  of  cloud 

SLY 

Left  end  of  cloud,  squared 

SLYS 

Right  end  of  Cloud 

SRY 

Right  end  of  cloud,  squared 

SRYS 

Hole  size 

SSH 

Hole  size,  squared 

SSHS 

B.6  SUBROUTINE  MPLACE. 

Subroutine  MPLACE  locates  each  burst  or  impact  point  of  the 
battery  volley.  The  aiming  and  precision  errors  are  computed  and  taken 
from  Subroutines  NORAN.  The  sequence  of  burst  placement  is  as  follows: 
For  volley  1  at  time  T-j,  weapon  1  fires  round  1  at  aimpoint  1;  weapon  2 
fires  round  2  at  aimpoint  2;  and  so  on  until  all  weapons  have  fired.  For 
volley  2  at  time  t2,  the  weapon  firing  is  repeated.  This  sequence 
continues  until  all  volleys  are  fired.  The  volley  number  uses  the  index 
I,  the  weapon  or  round  number  uses  the  index  J.  Each  of  the  NSAMP 
replications  called  in  the  MAIN  program  constitutes  a  different  firing 
occasion. 


Figure  12  shows  the  development  of  the  impact  points  dia- 
gramatically.  The  weapon  axes  can  be  displaced  from  the  reference  earth 
axes  by  an  angle  ANG.  In  the  weapon  axes  system,  the  target  is  located 
at  point  Xj,  Yj,  which  is  also  the  centroid  of  the  volley  aimpoints.  On 
one  occasion,  the  centroid  displacement  due  to  the  aiming  error  is 


42 


9-m 


oe 

//•nPAcr 


I 


I — 

I 

ywi/w  I 


f  I 
R  ' 


I 


-  -  -X  /XC.P.o\ 
^  V  tItO' 


pAoJSe-T'/J.£'(2i 

Ih^PAC-T 


XCA^XlOlS'KU  (a-)  \ 
;  (  yCA  +  Y  IDEA  l_  (1-3  ) 


s 


^  I  2. 

ymfAL(3)| 

«  i-i 


9 


/ 


3 


£>rs/«rp  cturr*. 

OP  /MPACr 


FIGURE  12  -  DEVELOPMENT  OF  AIMING,  PRECISION 
AND  SUBMUNITION  IMPACT  ERRORS 


43 


determined  by  Subroutine  NORAN  and  located  at  position  XCA,  YCA.  The 
impacts  of  each  of  the  volley  rounds  are  displaced  from  XCA,  YCA  by 
their  aimpoint  shifts  and  by  their  individual  precision  errors.  The 
latter  is  also  determined  by  Subroutine  NORAN.  The  final  locations  of 
each  round  is  XBA(l),  YBA(l);  XBA(2),  YBA(2);  etc.  Normally  ZBA  is 
zero.  The  time  each  volley  is  delivered  is  given  by 

BRSTTM  =  (I-l)  *  VT 

where  I  is  the  volley  number  and  VT  the  time  between  volleys. 

The  weapon  axes,  which  define  the  line  of  fire,  are  usually  selected 
coincident  with  the  earth  axes  in  which  case  ANG  is  zero  as  in  Figure 
12. 


B.7  SUBROUTINE  NORAN 

Subroutine  NORAN  determines  the  aiming  error  for  each  occasion 
or  program  replication  and  the  precision  error  for  each  round  delivered. 

B.8  SUBROUTINE  COMBST 

Subroutine  COMBST  transforms  the  impact  or  burst  points 
XBA(I),  YBA(I)  computed  in  subroutine  MPLACE  into  the  wind  coordinate 
system  with  an  origin  coincident  with  that  of  the  earth  fixed  system  and 
an  X-axis  in  the  wind  direction.  In  this  subroutine  the  index  I  covers 
all  rounds  fired,  WVUl  and  WVU2  are  the  direction  cosines  between  axes. 
The  choice  of  wind  axes  is  made  to  facilitate  the  integration  of  smoke 
densities  along  selected  lines  of  sight  in  Subroutine  ICON. 

B.9  SUBROUTINE  TIME 

Subroutine  TIME  is  executed  at  each  successive  period  MT, 
selected  in  MAIN,  for  a  total  of  NT  periods.  At  each  period,  the 
quantity  AGE  represents  the  age  of  all  volleys  or  bursts  in  the  program. 
The  function  MUNGRP  is  the  chronological  order  of  these  volleys  of 
bursts  and  is  used  as  an  index  for  identifying  successive  values  of 
burst  sigmas.  The  location  of  each  burst  at  period  MT  is  calculated 
in  the  earth  axes  system  and  given  as  XCENT(I),  YCENT(I),  where  the 
index  I  is  the  burst  number.  XCENTP(I)  is  the  location  of  the  burst 
point  in  the  wind  axes  system.  Finally,  the  maximum  component  cfXACH* 
°YA(n  of  each  burst  after  resolution  into  the  wind  axes  system 
is  calculated  for  each  program  time  interval,  MUNGRP. 
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B.IO  SUBROUTINE  SIZE 


Subroutine  SIZE  establishes  the  boundaries  of  a  rectangle,  in 
the  earth  axes  system,  which  encloses  all  bursts  to  be  delivered.  The 
rectangle  dimensions  are  determined  by  examining  the  ±4cjxa  and 
values  of  these  bursts  and  selecting  maximum  values.  This  rectangular 
area  serves  as  an  exploratory  region  for  determining  smoke  concentrations, 
{see  Figure  3). 

B.ll  SUBROUTINE  MXMIN 

Subroutine  MXMIN  determines  a  set  of  lateral  limits  y] ,  y2, 
which  are  used  for  establishing  the  array  of  observer  "look"  lines  in 
Subroutine  CALC.  These  "look"  lines  pass  through  the  ±3cjxa  coordinates 
of  the  extreme  bursts  of  the  cloud.  Since  the  slope  of  these  "look" 
lines  is  needed,  the  lateral  points  yj.,  y2,  are  calculated  at  a  selected 
distance  YTLINE  downrange  in  the  earth  axes  system. 

B.12  SUBROUTINE  MATCON 

Subroutine  MATCON  takes  the  ±4o  dimensions  of  the  smoke 
rectangle  from  Subroutine  SIZE  and  determines  a  40  X  40  array  of  x-y 
coordinates  within  these  limits.  These  coordinates  are  then  resolved 
into  the  wind  axes  system.  By  calling  Subroutine  DENSTY,  the  concen¬ 
tration  of  smoke  in  milligrams  per  cubic  meter  at  each  coordinate  is 
computed.  This  concentration  is  designated  by  the  symbol  C0NMAT(I,  J) 
which  is  the  same  as  CON  used  in  subroutine  DENSTY. 

B.13  SUBROUTINE  DENSTY 

Subroutine  DENSTY  determines  the  total  concentration  of  ob¬ 
scurant  in  milligrams  per  cubic  meter  due  to  all  bursts  at  each  of  the 
40  X  40  array  of  points  calculated  in  Subroutine  MATCOM.  The  Mass 
concentration  at  each  point  is  calculated  by  the  expression  given  in 
Equation  (2.3).  The  transport  of  the  cloud  is  determined  by  the  expression 

X  =  XP  -XCENTP 

where  XCENTP  is  the  burst's  centroid  location  measured  in  the  wind  axes 
system  and  obtained  from  Subroutine  TIME.  XP,  YP,  and  ZP  refer  to  the 
array  of  points  obtained  from  Subroutine  MATCON.  No  transport  is  de¬ 
veloped  in  the  YCENTP  and  ZCENTP  directions  which  are  normal  to  the  wind 
vel ocity. 

B.14  SUBROUTINE  MATCL 

Subroutine  MATCL  is  a  supplementary  procedure  for  examining 
the  concentration  of  smoke  agent  along  numerous  lines  of  sight  passing 
through  the  exploratory  rectangular  area  established  in  Subroutine  SIZE. 
Two  functions  are  performed  in  this  Subroutine,  they  are: 
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a.  Along  the  lower  boundary  of  the  rectangle,  40  points  are 
selected  each  of  which  will  be  used  as  the  source  of  a  family  of  "look" 
lines  extending  to  all  of  the  40  points  on  the  upper  boundary. 

b.  The  source  coordinates  of  each  line  and  the  line's  direc¬ 
tion  cosines  are  transformed  from  the  earth  axes  to  the  wind  axes. 

After  this.  Subroutine  ICON  is  called  to  determine  the  integrated  concen¬ 
tration  of  the  obscurant  along  each  line.  This  is  noted  as  CONLIN  (I,J) 
in  milligrams  per  square  meter. 

B.15  SUBROUTINE  ICON 


Subroutine  ICON  integrates  the  mass  of  obscurants  along  each 
"look"  line,  through  the  gaussian  density  distributions  of  each  burst 
for  each  time  period  MT.  The  "look"  lines  are  determined  in  Subroutine 
MATCL.  Mass  concentrations  along  specific  "look"  lines  originating  from 
a  fixed  observer  position  and  determined  in  Subroutine  CALC  also  employ 
ICON.  The  line  concentration  is  given  in  milligrams  per  square  meter 
and  is  computed  in  accordance  with  the  following  equation: 


NBURST 

'  I  .(0^(1)  0^1)  0^(1)  * 


X  exp 


1  r  /’x+S*Dx-xcent  ^ 


where  DX,  DY  and  DZ  are  direction  cosines  and  A  and  S  are  functions  of 
the  burst  centroid  location  and  its  standard  deviations  and  are  defined 
in  the  program  listing.  Wind  axes  are  used. 

B.16  SUBROUTINE  CALC 

Subroutine  CALC  determines  the  visibility  through  all  parts  of 
the  smoke  cloud  from  an  observer  location  XOBS,  YOBS.  Lines  of  sight 
are  selected  to  various  points  incrementally  spaced  a  distance  ZINC 
apart  and  located  at  a  distance  YTLINE  from  the  observer.  Along  each 
line,  the  mass  of  smoke  agent  is  computed  by  calling  Subroutine  LCON. 

Each  line  of  sight  which  permits  observation  through  the  cloud  contributes 
a  length  of  XINC  to  a  "hole"  size  in  that  cloud.  The  total  length  of 
the  cloud  is  defined  as  the  distance  between  the  extreme  point  through 
which  obscurance  first  occurs  and  last  occurs.  Figure  13  shows  the 
construction  of  the  observer-cloud  geometry  treated  in  Subroutine  CALC. 
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FIGURE  13  -  ANALYSIS  GEOMETRY  AS  USED  IN 
SUBROUTING  CALC 
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B.17  SUBROUTINE  VEVAL 


Subroutine  VEVAL  computes  the  left  and  right  extremes  of  the 
obscuring  screen  LY,  RY;  the  screen  length,  SLY;  and  the  accumulated 
hole  size,  SH  across  the  screen.  In  preparation  for  statistical  averaging, 
the  summation  of  the  square  of  these  quantities  is  determined  for  specific 
periods  MT  for  each  replication  of  the  program  KI.  The  initialization 
of  these  quantities  has  already  been  established  in  Subroutine  CLEAR. 

The  quantity  NS AMM,  which  defines  the  number  of  times  no  screen  is 
detected,  is  also  calculated. 

B.18  SUBROUTINE  CALPRT 

Subroutine  CALPRT  sums  up  the  main  results  of  the  program  by 
averaging  the  screen  length,  SSL;  "hole"  size,  SSH;  and  the  left  and 
right  coordinates  of  the  cloud,  SLY  and  SRY  respectively,  at  each 
specific  period  MT  for  a  total  of  KI  replications  of  the  program. 
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APPENDIX  C 


DESCRIPTION  OF  SUBROUTINES  FOR  HC  MUNITIONS 


C.l  GENERAL 

The  subroutine  for  HC  munitions  are  similar  to  those  for  WP 
munitions.  There  are  two  significant  differences.  First,  the  HC  smoke 
cloud  is  defined  as  a  series  of  small  periodic  bursts  or  puffs  origina¬ 
ting  at  the  point  of  the  munition's  impact.  The  number  of  bursts  and 
the  elapsed  time  between  them  is  defined  as  NPUFTL  and  TINCR  respec¬ 
tively.  Representative  values  for  these  quantities  used  in  the  past  are 
120  bursts  at  1  second  intervals.  Secondly,  for  those  cases  where  the 
munition  contains  a  number  of  submunitions,  the  latter  are  dispersed 
radially  in  accordance  with  a  selected  non-linear  distribution  and 
dispersed  directionally  according  to  a  uniform  distribution.  Flow  charts 
for  the  HC  program  are  given  in  Appendix  D. 

C.2  HC  PROGRAM,  MAIN 

MAIN  HC  contains  essentially  the  same  subroutines  as  MAIN  WP. 
These  similarities  include  provisions  for  calculating  the  concentration 
length  thresholds  for  various  sensors  (NDEVIC)  so  that  more  than  one 
device  can  be  handled  in  one  run.  Also  the  earth  fixed  axes  is  assumed 
to  be  aligned  along  the  line  of  fire  of  the  delivery  weapon  thereby 
eliminating  the  weapon  axes  system. 

C.3  SUBROUTINE  STCL 

This  subroutine  is  also  similar  to  that  of  the  WP  program 
except  that  the  fill  weight  per  puff,  rather  than  fill  weight  per  muni¬ 
tion,  is  considered.  Also  the  time  between  puffs,  TINCR,  is  identified. 

C.4  SUBROUTINE  CLEAR 

As  in  the  corresponding  WP  subroutine,  CLEAR  sets  initial 
values  of  statistical  variables  to  zero. 

C.5  SUBROUTINE  MPLACE 

In  addition  to  placing  the  munitions  in  accordance  with  firing 
accuracy,  the  HC  subroutine  also  places  the  submunitions  in  accordance 
with  the  radial  and  directional  distribution  of  these  quantities.  This 
subroutine  calls  NORAN  and  MUNDSP  to  identify  all  errors. 

C.6  SUBROUTINE  NORAN 

As  in  the  WP  program,  NORAN  determines  the  aiming  error  and 
precision  errors  of  each  round  delivered. 
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C.7  SUBROUTINE  MUNDSP 


Scatterable  submum’tions  contained  in  certain  projectiles  are 
placed  in  accordance  with  radial  and  directional  distribution  functions. 
As  already  stated.  Subroutine  MUNDSP  determines  these  locations  using  a 
nonlinear  and  uniform  distribution  respectively.  These  displacements 
are  used  in  Subroutine  MPLACE  along  with  the  aiming  and  precision  errors 
from  Subroutine  NORAN  to  define  the  impact  points  for  all  submunitions. 

C.8  SUBROUTINE  CONBST 

This  subroutine  is  identical  to  that  of  the  WP  program.  The 
impact  points  of  the  submunitions  are  transformed  from  the  earth  axes  to 
the  wind  axes  system. 

C.9  SUBROUTINE  TIME 

Subroutine  TIME  for  HC  munitions  has  been  expanded  over  the 
corresponding  Subroutine  for  WP  smoke.  As  before,  the  displacement  and 
diffusion  of  the  puffs  are  computed  as  a  function  of  time.  In  addition, 
the  location  and  size  of  the  first  and  last  puff  of  each  submunitions 
are  determined  in  order  to  define  the  extremeties  of  the  cloud  in 
Subroutine  SIZE. 

C.IO  SUBROUTINE  CON 

Subroutine  CON,  the  location  of  the  first  and  last  puff  from 
subroutine  TIME  are  transformed  to  the  earth  axes  system. 

C.ll  SUBROUTINE  SIZE 

Subroutine  SIZE  establishes  the  boundaries  of  an  exploratory 
rectangle  in  the  wind  axis  system.  The  rectangle  is  determined  by  the 
±4a  values  of  the  extreme  puffs  considering  all  of  the  submunitions  of 
the  delivered  load. 

C.12  SUBROUTINE  MXMIN 

As  in  the  corresponding  subroutine  for  WP  smoke,  this  program 
determines  the  lateral  limits,  now  called  XMIN  and  XMAX,  used  for 
establishing  an  array  of  observer  look  lines.  These  "look"  lines  are 
bounded  by  the  ±4a  coordinates  of  the  extreme  puffs  determined  from 
Subroutine  SIZE.  These  lateral  limits  are  calculated  using  their  slopes 
to  the  origin,  RATMN  and  RATMX. 

C.13  SUBROUTINE  MATCON 

This  Subroutine  is  identical  to  that  for  WP  smoke.  A  40  x  40 
array  of  x-y  coordinates  are  located  within  the  limits  of  the  exploratory 
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rectangle  taken  from  Subroutine  SIZE.  These  points  are  then  resolved 
into  the  wind  axes  system.  By  calling  Subroutine  DENSTY,  the  concen¬ 
tration  of  smoke  in  milligrams  per  cubic  meter  at  each  coordinate  is 
computed  and  noted  by  the  symbol  CONMAT. 

C.14  SUBROUTINE  DENSTY 

DENSTY  calculates  the  displacement  of  any  point,  XP,  YP,  ZP  in 
the  wind  axes  system  from  the  munition  burst  points.  DENSTY  then  calls 
Subroutine  DENM  to  determine  the  density  of  smoke  at  each  of  the  desired 
points  due  to  the  individual  effects  of  all  puffs  in  the  field. 

C.15  SUBROUTINE  DENM 

The  integration  of  the  individual  densities  of  all  puffs  in 
the  field  is  performed  by  Subroutine  DENM  at  the  points  of  interest. 

The  densities  of  all  puffs  at  all  points  are  computed  in  Subroutine 
DENPF. 

C.16  SUBROUTINE  DENPF 

The  density  of  smoke  at  any  point  x,  y,  z  for  each  puff  of 
each  submunition  is  calculated  by  this  subroutine. 

C.17  SUBROUTINE  MATCL 

This  subroutine  is  identical  to  the  corresponding  subroutine 
for  the  WP  program.  MATCL  is  a  supplementary  procedure  for  examining 
concentration  of  smoke  agent  along  numerous  lines  of  sight  passing 
through  the  exploratory  rectangular  area  established  in  Subroutine  SIZE. 
Two  functions  are  performed  in  this  Subroutine,  they  are: 

a.  Along  the  lower  boundary  of  the  rectangle,  40  points  are 
selected  each  of  which  will  be  used  as  the  source  of  a  family  of  "look" 
lines  extending  to  all  of  the  40  points  on  the  upper  boundary. 

b.  The  source  coordinates  of  each  line  and  the  line's  direc¬ 
tion  cosines  are  transformed  from  the  earth  axes  to  the  wind  axes.  After 
this.  Subroutine  LCON  is  called  to  determine  the  integrated  concentra¬ 
tion  of  the  obscurant  along  each  line.  This  is  noted  as  CONLIN  (I,J)  in 
milligrams  per  square  meter. 

C.18  SUBROUTINE  LCON 

LCON  computes  the  displacement  of  the  exploratory  point  x,  y, 
z  (called  POX,  POY,  POZ)  from  the  submunition  burst  point.  Using  this 
information  and  the  slope  of  the  LOS  from  the  observer  to  that  point, 
integration  of  all  smoke  densities  along  the  designated  line  is  per¬ 
formed  by  calling  Subroutine  LCONM.  The  concentration  along  the  LOS  is 
called  TOTLNC  in  these  Subroutines. 
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C.19  SUBROUTINE  LCONM 


Subroutine  LCONM  integrates  through  the  puffs  of  each  sub¬ 
munition  or  burst  along  the  LOS  defined  by  a  point  (POX,  POY,  POZ) 
and  the  slope  through  that  point  to  the  observers  position.  LCONM 
calls  Subroutine  LCONPF  for  supplementary  calculations  in  carrying  out 
the  integration. 

C.20  SUBROUTINE  LCONPF 

The  total  mass  of  smoke  along  the  observer's  LOS  is  determined 
by  integrating  through  the  Gaussian  trivariate  distribution  of  each  puff 
in  turn.  LCONPF  performs  supplementary  calculations  for  use  in  Subroutine 
LCONM. 

C.21  SUBROUTINE  CALC 

Subroutine  CALC  for  the  HC  program  is  identical  to  that  for  the 
WP  program.  CALC  determines  the  visibility  through  all  parts  of  the 
smoke  cloud  from  an  observer  location  XOBS,  YOBS.  Lines  of  sight  are 
selected  to  various  points  incrementally  spaced  a  distance  XINC  apart 
and  located  at  a  distance  YTLINE  from  the  observer.  Along  each  line, 
the  mass  of  smoke  agent  is  computed  by  calling  Subroutine  LCON.  Each 
line  of  sight  which  permits  observation  through  the  cloud  contributes  a 
length  XINC  to  a  "hole"  size  in  that  cloud.  The  total  length  of  the 
cloud  is  defined  as  the  distance  between  the  extreme  points  through 
which  obscurance  first  occurs  and  last  occurs.  Figure  13  shows  the 
construction  of  the  observer-cloud  geometry  treated  in  Subroutine  CALC. 

C.22  SUBROUTINE  ENDPTS 

ENDPTS  compiles  statistical  information  on  the  left  and  right 
extremes  (±4a  values)  of  the  exploratory  rectangle  defined  in  Sub¬ 
routine  SIZE.  This  information  is  used  in  Subroutine  CALPRT  to  deter¬ 
mine  the  average  location  and  standard  deviation  of  the  upward  and  down¬ 
wind  boundaries  of  this  rectangle  at  each  time  the  smoke  screen  is  ex¬ 
amined  (SNAP  time).  This  information  was  meant  for  future  use  of  the 
program  and  contributes  little  to  any  results  as  it  is  presently  listed. 

C.23  SUBROUTINE  VEVAL 

Subroutine  VEVAL  for  the  HC  program  is  identical  to  that  for 
the  WP  program.  This  subroutine  computes  the  left  and  right  extremes  of 
the  obscuring  screen  LY,  RY;  the  screen  length,  SLY;  and  the  accumulated 
hole  size,  SH  across  the  screen.  In  preparation  for  statistical  averaging, 
the  summation  of  the  square  of  these  quantities  is  determined  for  specific 
periods  MT  for  each  replication  of  the  program,  KI.  The  initial  quantities 
have  already  been  established  in  Subroutine  CLEAR.  The  quantity  NSAMM, 
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which  defines  the  number  of  times  no  screen  is  detected,  is  also  calcu¬ 
lated  and  used  in  adjusting  statistical  computations  in  Subroutine  CALPRT. 

C.24  SUBROUTINE  CALPRT 

The  main  results  of  the  program  are  given  in  CALPRT.  These 
are  obtained  by  averaging  the  screen  length,  SSL;  "hole"  size,  SSH;  and 
the  left  and  right  coordinates  of  the  cloud,  SLY  and  SRY  respectively, 
at  each  specific  period  MT  for  a  total  of  KI  replications  of  the  pro¬ 
gram.  In  addition,  CALPRT  uses  information  from  Subroutine  ENDPTS  to 
determine  the  probablelocations  of  the  exploratory  rectangle  established 
in  Subroutine  SIZE. 
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vr 
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V/AP 
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//JAMP 

vr 
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CALL  STCL  (S^  0^,(2 
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lAJVdL  ^VdT 
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curoFU  ^-YrL»Kyr 
coroFi®^  ■*•  yTLihiff 
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S^O'3  ^  O  T’/a^  ^ 
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Sc'a  O  (^>r.  t>r,  ^yoL,  a/s  } 
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SucpouriMiT  MXMINJ  (  D,'i\  .yz-) 


66 
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67 


Sucir3^or/Ay£’  OEM  STY  (  xp,  yp,  iPy  co/^) 


68 


ATC.1^  xhi^Mj  vlow, 
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So>S£our-//og-  dA^LC.  ( xo3S,  Voss^  ynif-'S’^  xmim,  xfiAXj  cith/is) 
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Su3bout/kJ^  LCO/J  (  !  Poy  ^  ,  oy , 


SoBi&oor/Aj£‘  VS'VAL  ^I-O 
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MAIN  HC 
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vr^  SWAP 
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Call  STCL  (qwukJ  .yiELO,  Bwa-tor,  3,  u) 


Do  is*i  1*1,  porvic. 


CLTHRS(I5  =  ALO^^|-TH£e'SCl.))/<l-CArTNtl.)J 
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=  AJVOL*  4  /s/mPA. 


o 
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Do  IT^  I,  NT 


T  ‘  Float  (IT)  *  SNAP 


Call  TIME  (:t) 

Call  Com 

Call  Sl^E  (^xlow,  xhibh^  vLovt/^ymovn) 
Call  mxmim  (  vthns^  xmim,xi^ax) 
Call  MAT com  xLow^xHi&H,yio\./,yHic,H) 

Call  MATCL  (^2,xLovv',.xHiGW^>fci3s^ysLiN£) 

1  rAi  r  /'Xoas.yoBi,  »,  yn  iwf,  xmax 
C-ALL  L-/M-L- ,C.LrH^S  ) 

Call  EMDPT5^xMi»j^xMAXjrT) 


I*  Do  2.‘i‘?  ze.  -  I ,  NOBVi  c 


Call  VEVA1_  (it,  tr.) 


Call  clear  Cdl.) 

Do  looo  K1 

-  Ij  MiAMP- 

^  r  NSAM^.  A/r,SNAP 

Call  CAlPRT^  A/os^tc-  j 

STO  p 

eKJ  O 
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Sunaour INS  M PLACE  (^nvoLj  nn^pa^  yr,  yr^yr 
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98 


Susnoun  J-C  ^coNn^uKo) 
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SoSftoOTloe  LCOh^PF  C  Mut^P  ,  tP  ,  pFLhJC) 


SufiffcouTiMe  Calc  (xoas.  yoqs,  so»vic,clthrs) 


CALL  LCOM  ^xoas  , 

.  DXP,  Dyp,  D»,  ToruMC.  ) 


DX  =  Ox/VI^A(i 
DY  =  0Y/VMA6,. 
Ol  =  D^/V/^A<5 


DXP  =  wYL/i  •*  DX  +  wvui  «  oy 
DVP  *  -  w  vD  i  ■»  DX  +  WYUI  ■»  OY 
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^uti/ijooT“i £NDPT3  ( , /r) 
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SuanourtNe  \/£VAL 
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Suaacor ffjg  ^^VAL  (cat^r.^ 
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SuQfiouTi^s  '^^VAL  (cot^r.) 


Suao.ouriA>e  C.ALP/ZT  ^ajap,  ^vou^vr^/^ofv/a) 
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Su3fZQuTif^e  CALPRT.  (cor^T.) 
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PAGE 


I 


PROGRAM  SMOKEWPC  INPUT, OUTPUT, TAPEH-INPUT) 

DIMENSION  TITLECBI .CATTN(3) .THRES(3) ,CL  THRS(3) 

COMMON  fP£Lf  X0(6dJ  , YB(60) ,ZB(BOI ,T0(6O ) ,NBURST 
COMMON  /RELB/XCENTPtBO) ,YCENTP(BOI 

COMMON  /WPSCR/  SIOX (50) ,SI0y(50) .SI0Z(50) 

COMMON  /WPSCR/  FACT T( 50 1 

COMMON  /WPCLD/  FSIOX,FSIOY,FSIGZ 
COMMON  /WPCLD/  DSIOX,DSIOY,DSIGZ 
COMMON  /WPCLD/  ASIOX,ASIOY,ASIOZ 
COMMON  /V^’CLD/  ESI0X,ESI0Y,ESI0Z 
COMMON  /WPCLD/  WNOSPO, FACTOR 

COfiON  /ABS/XBA(60) ,YBA(BO) ,ZBA(BO) .TBA (60) ,NBRSrA 
COMMON  /ABSe/  XCENT (BO) .YCENTIBO) ,ZCENT (60) 

COMMON  /ABS3/  tlGXA (50) .SI0YA(50) 

COMMON  /TMSCLE/INSNAP, INVOL.DT 

COMMON  /PLACE/  SI06R,SI080,SIOAR,SIGAD. REL 

COMMON  /PLACE/  XIDE AL ( 6) , YIDEAL (6) ,ZIOE AL (6) 

COMMON  /kfODIR/WVUI  ,WVUB 
COMMON  /AGCINO/  UOE(BO) 

COMMON  /FLDSIT/  CUTOFL.CUTOFR.XINC 

DATA  FSIGX/ . 1553/ .FSI0Y/3.m / ,FSI OZ/ I .35/ 

DATA  DSIOX/I .0/ ,DSI GY/ 100. 0/ ,DSI0Z/30. 0 / 

DATA  ZCENT/50»0.0/ 

500  CONTINUE 

READ(‘4,3)  TITLE 
3  FORMAT! BAIO) 

PRINT  10, TITLE 
10  FORMAT! IHI .BAIO) 

READ(H,  t  )U, ALPHA, BETA 
I  FORMAT (3r 1 0.5) 

PRINT  ‘♦.U,  ALPHA, BET  A 

H  FORMATdHO,*  WIND  SPEEDtX  DIRECT)  -*,FI0.5, 

•  ‘ALPHA  -•.FI0.5.*BETA  -*,FI0.5) 

READ(H,II)  UVUI.UVU3 
n  FORMAT(3FI0.5) 

PRINT  I3,WVUI.WVU3 

13  FORMAT! I  HO, ‘HI  NO  VECTOR  DIRECTION  COSINES  X  -*,Ft0.5.‘Y  -‘.FIC.B) 

READ(H,6)  YIELD,QMUN,SIGXS,SIOYS.SIGZS 

6  F0RMATtFI0.5,EI0.5, 3ri0.5) 

PRINT  3.YIELD,QMUN,  SIGXS,SIDYS,SIGZS 
3  F0RMAT(*0*,‘  YIELD  -* ,FI 0. 5. *OUAN  OF  MUN  (ADJ)  -‘.EIC.B. 

I  *  SOURCE  SIGMAS  •‘.3FI0.5) 

READ(H,5)  NOEVIC 
5  FORMAT! 15) 

READ(H,7)  ( (CATTN( I R) ,THRES( IR) ) , I R- 1  ,  NDEVI C ) 

7  FORMAT (BF 10. 5) 

CALL  STCLIQMUN. YIELD. SIGXS,SIGYS.SIGZS. ALPHA. BETA. U) 

DO  159  !•! .NCCVIC 

CLTHRS! I )-ALOO( I .0“THRES( I ) )/ (-CATTN( I ) ) 

PRINT  95. I .CATTN( I ) .THRES( I ) ,CLTHRS( I ) 

95  FORMATCO‘.I3.‘  COEF  ATTN  -‘.FI0.5. 

I  *  THRESHOLD  -‘.F10.3. 


Ill  A 


PAGE  1 A 


?  •  CL  THRES 

159  CONTINUE 

C  . . MUNITIONS  PLACEMENT  It  RE 


READti+.SO)  NRPV 
50  FORMAT!  MO) 

PRINT  m.NRPV 

!»♦  FORMATdHO, •NUMBER  OF  ROUNDS  PER  VOLLEY  *,14/ 

2  IHO,* IDEAL  IMPACT  POINTS  •/ 

•  *0*  .  •  NO  .  •  ,6X,  'X*,  MX,  'Y*  .  1 1  X,  -ZM 

DO  2500  J-l .NRPV 

REAOtU.SS)  XIDEAL(J) , Y IDEAL fj) .ZIDEALIJ) 

55  FORMAT! 3F I  0.5) 

PRINT  16. J.XIDEAL! J) .YIDEAL! J) .ZIDEALtJ) 

16  FORMAT!*  •,m.6(IX.F10.2.1X)) 

2500  CONTINUE 

R£ADtM.2M  NSAMP.NVOL.NT.OT.X 

21  FORMAT(3I5.2F5.0) 

PRINT  23,NSAMP,NVX  .NT.DT.X 
23  FORMAT! IHO, ‘TOTAL  NUMBER  OF  SAMPLES  -*,15./ 

1  IX. ‘TOTAL  NUMBER  OF  VXLEYS  «‘,I5,/ 

2  IX. ‘TOTAL  NUMBER  OF  TIMES  15./ 

3  IX. ‘Tire  INCREhCNTS  -•  ,F5. 0, ‘SEC  .  ‘  ,  / 

IX. ‘LINE  INCREMENTS  *‘,F5,0) 

REAOm.lB)  VT.SNAP 

18  FORMAT12F10.  I  ) 

PRINT  I  9. VT.SNAP 

19  FORMAT!  IHO, ‘TlhC  BeT^€EN  VXLEYS  -‘.FlO.l./ 

1  IX. ‘SNAP  Tire  INCRErCNT  -‘.FIO.M 

READC+.aa)  SIGBR.SIGBO.SIOAR.SIGAD.REL 

22  F0RMAT(5F10.4> 

PRINT  e4,SIG8R.SIG0D.SlGAR.SlGAD.REL 
2H  FORMAT! IHO. ‘BR  - ‘ . F 5 . 0 . 5X , ‘BO  -‘.FS.O,/ 

1  IX. ‘AR  -‘.F5.0.5X.*AD  -‘.FS.O./ 

2  1X.‘REL  -‘.FS.O) 

NBRSTA«NVX*NRPV 

tNSNAP-SNAP/DT 

[NVX-VT/OT 

CALL  FRMCLDfNT.DT.NVOL.NS) 

XOBS-0.0 

''OBS-0.0 

READ(4.20)  YTLirC.CUTDFL.CUTOFR 

20  FDRMAT!3FI0.2) 

REAO!‘4.20)  XT.YT 
READIH.aO)  XINC 
Z-l  .5 

PRINT  73.X08S,Y08S 

73  FORMAT! IHO. ‘OBSERVER  COORD  -‘.BFIO.M 
PRINT  71 .XT.YT 

71  FORMAT! IHO. ‘AIMPOINT  COORO  -‘.BFIO.M 
PRINT  72.YTL1^C 

72  FORMAT!  I  HO. ‘DISTANCE  TO  LirC  •‘.FIO.M 
PRINT  7»4.XINC 

7H  FORMATMHO,*  X-INCRErCNT  -‘.F5.M 
PRINT  75.CUT0FL.CUT0FR 

75  FORMATMHO.*  90  DEGREE  SECTOR  -•.F7.0.*  .  •.F7.0) 
PRINT  35. Z 

35  F0RMAT!*0*. ‘HEIGHT  ABOVE  CENTROID  ‘.FI0.5) 

CALL  CLEAR  fX) 

DO  28  K 1 « I . NSAMP 
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PRINT  8,KI 

B  FORMAT  I  MV  DATA  FOR  SAMPLE  NUMBtP  Ml**! 

CALL  MPLACECNVOL  ,NRPV.VT,XT,YTJ 
CALL  CONBST 

PRINT 

HP  FORMATC-OM*  MUNITION  BURSTS  AT  M 
PRINT  H3 

H3  FORMAT  CO-  .  I  IX.  -LOCATION*  ,  I  I  X ,  5X ,  7X  .  •  T  I  ME  •  ) 

I-O 

DO  199  IV*1 ,NVOL 
DO  199  IRD-I .NRPV 
I-l 

PRINT  HB. I . I V, IRD.XBAC  n , YBA( I  1 .ZBAI I J . TBA( I ) 

HB  FORMATCO-.  3(IX.I3.IXJ  .3(IX.F8.P.IX».6X.F8.P» 

199  CONTINUE 

DO  I  100  MT-I .NT 

CALL  TIMEIMT.UJ 

C  PRINT  3P 

3P  FORMATI-0-.-  MUNITION  CHARACTERISTICS  •> 

C  PRINT  33 

33  FORMAT I-O- . I  I X. -LOC  AT  ION- . II X .5X , IHX. -S I ZE- > 

DO  350  IMUN-I .NBRSTA 
MG-I AOEI IMUN) 

IF(MO.LE.O)  GO  TO  360 

C  PRINT  38.  IMUN.XCENT  I  IMUNl  .YCENTI  IMUN)  .YCENTI  IMUN)  . 

C  *  SIOXAI  IMUN)  .SIOYAI  IMUN)  .S)CZ(MO) 

38  F0RMAT(-0-. IX. IP.PX .31 IX.FB.P. IX) .5X. 31  I X .F8.P . I  X) ) 

350  CONTINUE 
360  CONTINUE 

CALL  SIZEIXLOW.XHIGH.YLOW.YHIGH) 

PRINT  86.XL0W.XHI0H  .YLOW.YHIGM 
86  FORMATIIHO.-  SCREEN  DIICNSIONS  X-- .PF7 .P .PX . - Y-- .PF7 .8) 

CALL  MXMINIYTLIKE.XniN.XMAX) 

PRINT  9.XMIN.XMAX 

9  FORMATIIH  .-XMIN  •- .F 1 0 .P.PX . -XMAX  --.FIO.P) 

YSLINE-500.0 


C  CALL  MATCONIZ.XLOW,  XH I GH.  YLOW.YHIGM) 

C  CALL  MATCL(Z.XLOW.XHIGH.YOBS.YSLINE) 

CALL  CALCIXOBS.YOBS.Z.YTLIKE.XMIN.XMAX.NOEVIC.CLTHRS) 
CALL  ENDPTSIXMIN.XMAX.MT) 

DO  899  IR-I .NOEVIC 
CALL  VEVAL  (MT.IR) 

P99  CONTINUE 
HOP  CONTINUE 
PB  CONTINUE 

CALL  CALPRT  (NSAMP, NT .SNAP.NVOL . VT .NOEV 1C ) 

STOP 

END 
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SUBROUT  I  CUEAR(a-) 

COMMON  /STATS/  SSL ( 50 . 3J .SSLS(50, 3) . SL Y J 50, 3> ,SL YS( 50, 3) , 

I  SRY(50,3» ,SRYS(50, 3» ,5SH(50,3J ,SSMS(50 .3J ,NSAMMf50.3) ,LENf  100) . 
a  O5L(50,3.I00) 

COMMON  /STAT?/  EXTXL ( 50) .EXTXLSC 50 ) ,EXT XR( 50 ) ,EXTXRS( 50 ) . NI  X ( 50 ) 

00  800  I -I .100 
800  LEN(  I  )«IX*  [ 

00  lao  I-l ,50 
NIX(  D-O 
EXTXLCn-O.O 
EXTXLS(I)-0.0 
EXTXRI I) -0.0 

EXTXRSt  n-0.0 
00  no  IR-1,3 
NSAMH(I.IR)-0 
SSL(I.IR)-0.0 
SSLS(I,IR)-0.0 
SLY(I,IR)-0.0 
SLYS(I,IR)»0.0 
SRY(I.tR)-0,0 
SRYS(I.IR)-0.0 
00  105  K- 1 . 1 00 
105  DSL( I . IR,K)-0.0 
SSH(I.IR)-0.0 
HD  SSHS(I,IR)-0.0 
rao  CONTINUE 
RETURN 
END 
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SUBROUTIfC  MPLACEtNVOL,NRPV,VT,XT,Yn 

COMMON  /ABS/  xereo.Yeieoi, 20(60), rereoi.NBURST 
COMMON  /PLACE/  SIGBR.SIGBO.SIOAR,SIGAD.  REL 
COMMON  /PLACE/  XIDE  AL  (6)  ,  YIOCAL  (61  .ZIDE  AL  (6) 

CALL  NORAN(R,5IOAR,E,SIOAD1 

XCA-XT+E 

YCA-YT^R 

DO  32  I-l.NVOL 

BRSTTM-  (I-M-VT 

DO  30  J-I,r«PV 

ir(RANF(DUM» .OT.REL  I  00  TO  30 
CALL  N0RAN(R,5I0eR,E, 510001 
XCRO-XCA^XIOEAL(J14E 
YCRO- YCA+ Y IDE AL ( J» ♦ R 

NUNNO-(  I-M-NRPV  J 

XB(MUNN01-XCR0 
YB(MUNNO>"YCRD 
2B(MUNN01"0.0 
TB(MUNN0I -BRSTTM 
30  CONTINUE 
32  CONTINUE 
RETURN 
END 
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MNQUTItC  CONBST 

CC3MM0N  /RCL/  XBC60)  ,  Y8C60)  ,Zec 60) ,  TB(60  )  .NBUKST 
COMMtJN  /ICL3/XCENTPe60)  ,YCENTP)60) 

COmON  /AB6/XBA(6QI  .YBA(60)  .ZBA(60)  .  TBA  C60)  .NBP6TA 
COmON  ./IMODIB/HVUI  ,WVU6 

nbmt-nbbsta 

DO  100  l-I.NBPSTA 

XB(  n«MVUt  •XBAt  I )  ♦NVUe*YBA(  t  I 

YB C  M  — -wvue* XBAe  M  ♦  WVUI  •  YBA ( I ) 

ZBtn«ZBA(M 
TB(  n«TBAt  I ) 

YCENTPCM-YBCI) 

100  CONTINUE 
MCTUNN 
END 
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SU0ROUTI^C  NORAN(R.  SR.O.SO) 
I  X-RANF(DUM) 

inx.LE.o.oi  GO  ro  i 

A«SQRTI-?.0'AL0GCX1 1 

e»6.ae3te53O7i0'RANr(ouM) 

R’A«SR*SIN(ei 

O-A*SD*COS(0J 

RETURN 

END 
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SUBROUTINE  S17ECXIL  .XUL.YLL,YUL) 

COMMON  /AB5/XBA(60) .Y0A(6O) .ZBA(60) .IBA (601 .NBRSTA 
COMMON  / ABS5/  XCENT ( 60 ) . YCENT ( 60  > . ZCENT  1 60 ) 

COMMON  /ABS3/  SI GXA ( 50) .SIGYA(50) 

COMMON  /AGEING/  I  AGE (60) 

XLL-  I .E6 
XIJL--1  .E6 
YLL-  I.E6 
YUL— I  .E6 

DO  BB  MUN-I .NBRSTA 
MUNGRP-IAOE(MUN) 

IF(MUNGRP.LE.O)  GO  TO  99 
TEST -XCENT ( MUN) . 0 ‘S I GXA ( MUN ) 

IF(TEST.LT.XLL)  xll-test 
TEST-XCENT  (MUN)  .  0  •  5 1  GXA  ( MUN ) 

IFdEST.OT.XVJL)  XUE-TEST 
test- YCENT  (MUN)  'H .  0  •S I OYA(  MUN) 

IF( TEST. LT. YLL)  YLL -TEST 
TEST-YCENT(MUN)+H.O*SIOYA(MUN) 

IF ( TEST. OT. YUL)  YUL -"EST 
B9  CONTINUE 
99  CONTINUE 
RETURN 
END 
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SUBROUTINE  MXMI N( YT L I NE ,XMI N. XMAX » 

COMMON  /WPSCR/  SIGX (501 .5I0YC50) .SICZtSOI 
COMMON  /ABS/X0A(6OI . YBA(60» ,ZBA(60) , TBA (601 .NBRSTA 
COMMON  /ABSa/  XCENT (601 .YCENT(60) .ZCENT (601 
COMMON  /AOeiND/  IAGE(60) 

RATMX--I .OE+IO 
RATrt4-I  .OE  +  IO 
DO  I  I-I,NBRSTA 
M-IAOCd) 

iriM.LE.O)  GO  TO  a 
ir(YCENT(n  .EQ.O.OI  GO  TO  I 
RArL-(XCrNT(  I  *5  IGX(MI  )/YCENT(  I  I 
ir(RATL.LT.RATMNI  RATMN^RATL 
RATN-( XCENT ( I l♦^. *5IGX(M) ) /YCENT( I  I 
ir(RATN.OT.RATMX)  RATMX-RATN 
t  CONTINUE 
a  CONTINUE 

XMIN-YTLI(C*RATMN 

XMAX-YTLINE*RATMX 

RETURN 

END 
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SUBROUT  I^r  CALC(XOBS,YQBS,Z,YTLINr  .XMlN.X^Ai.NpeVIC.CLTHRSl 

COMMON  /HNDOIR/WVUI  .WVUe 

COMMON  /FLDSIT/  CUTOf L .CUTOf R.XINC 

COMMON  /LINES/  V(  I  001 , 31 , I XLEFT , IXRGHT 

DIMENSION  CLTHRSrSl 

XL-AMAXICCUTOFL.XMINI 
XR-AMINI (CUTOFR.XMAX) 

IXLEFT-INT( (XL-CUT0FL)/XINC1  ♦  I 
1XRGHT»INT( (XR-CUTDFLl/XINC)  ♦  I 
DO  31  IP-IXLEFT.IXROHT 
X-XINC  •  FLOATCIP-I)  ♦  CUTOFL 

Dx*x-xoes 

DY«»YTL1NE-Y0B5 

DZ-0.0 

VMAG-SQRT (DX« ‘a^DY* 'a+DZ* -ai 

DX-DX/VMAG 

DY-DY/VMAO 

DZ-DZ/VMAO 

DXP-WVUI •Dx>wvua*DY 

DYP--WVUa«DX*WVUI 'DY 

CALL  LCON (XOBS. YOBS. Z.DXP.DYP.DZ, TOIL NC  1 
DO  a9  IR-I .NDEVIC 
Vt IP, IRJ-O.O 

IFdDTLNC.OE.CLTHRS  t  IRII  VC  IP,  IR)-I  .0 

ag  continue 

31  CONTINUE 

DO  35  IR-I .NDEVIC 

C  PRINT  3H.(V(1.1R1.1-IXLEFT,IXR0HT1 

3M  FORMAT C« 0*.  CaSFH.  I  . /n 

35  CONTINUE 
RETURN 
END 
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SL^RCXJTIhC  VEVALdf.lRl 


COMMON  ^FLDSIT/  CUTOTL.CUTOTR.XINC 
COMMON  /LINES/  V( 100 1 , 31 , I XLEFT , I XRGHI 

COMMON  /STATS/  SSL ( 50 . 31 ,SSLS(50 , 31 . SLY (50 . 3) .SLYSf 50 , 31 , 

I  SRY(50.3) .SRYS(50. 3) .SSH(50.3) ,SSHS(50 ,31 ,NSAMM(50.31 .LEN(  100)  . 
a  nsL(5p.3. 1001 

INTEGER  RY,SH 

LY»0 

RY-0 

SH-0 

SL-0.0 

C  LEFT  Y-COORD  •=  LY 

IF( IXRGHT'IXLEFT.LE .0)  GO  TO  95 
DO  10  I'IXLEFT, IXRGHT 
IFrV( I , IR) .EQ.O.Ol  GO  TO  10 
LY*XINC*FLOAT( l-I 1  ♦  CUTOFL 
K-I 

GO  TO  15 
10  CONTINUE 
GO  TO  MO 

C  RIGHT  Y-COORD  -  RY 

15  DO  aO  J-K, IXRGHT 

IF(V( J. IRl .EO.O.Ol  GO  TO  aO 
RY-XINC-FLOAT( J-I  1  ♦  CUTOFL 
M-J 

aO  CONTINUE 

C  TOTAL  SCREEN  LENGTH  -  SL 

SL-RY-LY 

IF(SL.EO.O. )  GO  TO  MO 
SSL ( M . IRl -SSL  ( I  I  , I R1 ♦SL 
SSLS( M . IRl-SSLSt 1 1  .IR)+SL«SL 
SLY(M.IR)-SLY(II.IR1+LY 
SLYS( I  I , IR1-SLYS(  I  I  , IR)+LY«LY 
SRY( I  I , IR)-SRY( M ,  I  R) ♦RY 
SRYSm  ,  IRl-SRYSt  M  .IR)^RY*RY 

C  INTERIOR  HOLE  LENGTH  •  SH 
C 

DO  30  I*K.M 

IF(V( I , IRl .EO. I .01  GO  TO  30 
SH-SH^XINC 
30  CONTINUE 

SSH( I  I , IR1-SSH(  I  I  .  I  Rl^SH 
SSHSt  I  I  .  IR1-SSHS(  I  I  ,  IRdSH'SH 
GO  TO  50 

C  NSAMM  TO  BE  USEO  TO  MOOIFY  NSAMP  SUCH  THAT 

C  AVE  LY.RY.AhO  SH  WI LL  BE  CALC.  GIVEN 
C  NON-ZERO  SCREEN  LENGTH. 

MO  NSArti(  I  I  ,  IR)-NSAMMt  I  I  .  IRdI 

C  DISTRIBUTION  OF  SCREEN  LENGTH 

50  DO  60  I -I , 100 

IF(SL.Ge:.LEN(  I)  )  GO  TO  60 
OSL(  I  I  .  IR.  I  )-OSL(  I  I  .  IR.  Idl  .0 
GO  TO  70 
60  CONTINUE 
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70  CONTINUF 
Rf TURN 
95  CONTINUE 

NSAMMt I  I . IRi=NSAMM(  I  I . IR»  ^  I 

RETURN 

END 
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SUBROUTINT  ENOPTS(XMIN.XMAX, IT) 

COMMON  /STATa/  EXTXL *50 ) .CXTXLSt W) » .tXT XRC 50 ) .EXTXRSr 50 ) . Nl  X(50 ) 
IFCXMAX-XMIN.LE.J.O)  00  TO  P5 
EXTXLC IT)-EXTXL» IT)  ♦  XMIN 
EXTXLSC IT)-EXTXLS(  I  T)  ♦  XMIN»XMIN 
EXTXRI|T)-EXTXR(  IT)  ♦  XHAX 
EXTXRS* IT)-EXTXRS( I T)  ♦  XMAX«XMAX 
RETURN 
P5  CONTINLC 

NIXC|T)-NIXeiT)  ♦  I 

RETURN 

END 
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SUBROUTINE  MATCONIZ  ,Xl  OW.XHIGH,  YL  OW.  THl 

COMMON  /WNDOIR/WVUl .WVU? 

OIICNSION  XCOROC^tO)  .YCORQfHOI 
01^CNS!0N  COMMA T (1*0  .‘♦01 


DO  1100  t-I ,H0 

XCOROt  I  I-ELOAT  fIM 
MOO  CONTINUE 

DO  1500  !«1  .‘♦0 
YcoRoin-rLOATM-n 


(I  .O/i+O.OMIXHIGH-XLOWMXLOW 


•(I  .O/i+O.OMlYHIOH-YLOWI+YLOW 


1500  CONTINUE 


PRINT  22 

22  FORMAT (•  X  *.5X,20X,».  T 

PRINT  26. (YCORDt IPi . IP-I .201 
26  rORMAT(»0»,7X,20(r6  .  I )  ) 

DO  2H00  I- I  .HO 
DO  2H00  J-I .HO 

X-WVUI  •XCOROI  H+WVU2»YC0R0t  J1 
Y-  -WVU2»XC0R0(  I )  ♦WVUI  -YCORDl  J1 
CALL  DENSTYIX.Y.Z.CONl 
CONMATI I , J)-CON 
2H00  CONTINUE 

DO  2500  1*1  .‘♦0 

PRINT  2e,XC0R0( 1 1 . ‘ CONMATI I , IPl , IP-I ,20 l 
IBLAM<-5»M/5)-I 
IF  (IBLANK.EQ.OIPRINT  32 
2500  CONTINUE 
PRINT  22 

PRINT  26, (YCOROt IPl . IP-21 -HOI 
DO  3000  I-I .HO 

PRINT  28.XCOROt  I  1  ,  tCONMATt I . IPl , IP-21 ,H0l 
IBLAM<-5»M/51-I 
IF  t IBLANK.EOOIPRINT  32 
3000  CONTINUE 

28  FORMAT!*  *, IX,F6. I .20tr6.0) 1 

32  FORMAT (•  *1 
RETURN 
END 
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SUBROUTINE  MATCL  t  Z  ,  XLOW.XHIGH,  YLOW,  YHIGIG 

COMMON  /WNDOIR/WVU!  ,WVUe 
DlfCNSION  XI tHO) .xaiHO) 

DlfCNSION  CONLINIMO .HO) 


YI-YLOW 

Ya-YHIGH 

DZ*0.0 

00  I  100  1-1  .HO 

XI  t  I  ) -FLOAT  (  I  -  I  )M  I  .O/HO.O)  MXHIOH-XLOW  1+XLOW 
I  100  CONTINUE 

DO  1500  I-l ,H0 

Xai I ) -FLOAT ( I-l ) M  1  .O/HO.O)* (XHIGH-XLOW) ♦XLOW 
1500  CONTINUE 


PRINT  I5,YI,Ya 

15  FORMATr»0«.«  Yl  -.FB.a.*  YB  '.FO.?) 

PRINT  aa 

aa  FORMAT (•  XI  •,5x,aox,»  xa  •) 

DO  8500  l-I .HO 
DO  8H00  J-l .HO 

DIST-SQRT(  (Xar J)-X 1 ( I ) ) ••8+rY8-Yl ) ••?) 

DX-(X8rj)-Xni))/DIST 

DY-IYe-YD/DlST 

X-WVUMXI  (  n  ♦WVU8»  YLDW 

Y— wvua»xi  (  I )  ♦WVUI  •  YLOW 

DXP-WVUI •DX+WVU8»DY 

DYP--WVUe«DX+WVUI -DY 

CALL  LCONIX.Y.Z.DXP.DYP.DZ.CON) 

CONLINI I.J)-CON 
8H00  CONTINUE 
8500  CONTINUE 

PRINT  86. rX8tI) , 1-1 .80) 

DO  89  K-|,H0 

PRINT  88. XI (K) . tCONLlNIK. I ) , 1-1 .80) 

89  CONTINUE 

PRINT  I5.Y1,Y8 
PRINT  aa 

PRINT  86.rx?(|),I-?|.H0) 

DO  87  K-1 .HO 

PRINT  aS.XI (K) . rCONLINlK, I) . 1-81 .HO) 

87  CONTINUE 

86  FORMAT (• 0« .7X. 801 F6 . 1 )  ) 

88  FORMAT!-  • ,F6. I ,80! F6.0) ) 

RETURN 

END 
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SUBROUTirC  CALPRT(NSAMP,NT.SNAP.NVni  .vl  .NDEVIC) 
r  CALCULATE  ANO  PRINT  STATISTICS 

COMMON  /STATS/  SSL ( 50 . 3 1 , SSLS t 50 . 3 ) . SL Y f 50, 3 ) , SL YS f 50 , 1 ) . 

1  SRY(50.3) .SRYSf50. 3) ,SSHt50.3) ,SSHS(50  ,31 ,NSAMMt50,3i ,LFN(  (00)  , 

2  OSL(50.3, 100) 

COMMOt^  /STAT^/  EXTXL(50) ,FXTXLS(50) ,EXTXR(50) .FXTXRS»50) ,N( X(50) 
PRINT  510.NSAMP 

510  FORMAT! IHI , •TOTAL  NUMBER  OF  SAMPLES  -‘.IS) 

DO  505  I -I .NT 
DO  520  IR«I .NDEVIC 
NSAMPI-NSAMP-NSAMM(  I  ,  IR) 

IFCNSAMPI  .LE.  ntX)  TO  519 

SLYS(  I  .  IR)-SQRT(  (SL  YS(  I  ,  1R)-SLY(  I  ,.IR)  “a.  /NSAMPl  )  /(NSAMPI  -  I  )  ) 

SRYSC  I  .  IR)-SQRT(  (SRYSt  I  .  IR)~SRY(1  ,  I R)  •  •  2 . /NSAMPI)  /  ( NSAMPI  - 1  )  ) 

SSHSI I  ,IR)-SQRT((SSHS(  I  .  IR) -SSHt  I ,  I R)  “a . /NSAMPI  )/( NSAMPI  -  I  )  ) 

SLY( I . IR)«SLY(  I . IR) /NSAMPI 
SRY(  r.  IR)-5RY( I .  IR) /NSAMPI 
S5H( I . IR)-SSH( I  ,  IR) /NSAMPI 

515  SSLS! I, IR) -SORT!  (SSLS( I , IR)-SSL< I , IR) •  •  2 . /NSAMP) / ( NSAMP- 1 ) ) 

SSL(I . IR)-SSL( I .  IR) /NSAMP 
520  CONTINUE 

NIXI-NSAMP-NIX(I) 

IF(NIXI .LE. I )  00  TO  505 

EXTXLSC  I  )-SQRT(  (EXTXLSC  I  l-EXTXLi!  I)**2./N[XI)/fNIX)-n) 

EXTXRSI  I  )-SQRT(  (EXTXRSC  I  )  -EXTXRi!  I  )  •  *2,  /  N I X I  )  /  (NI X I  -  M  ) 

EXTXLI  n-EXTXL(  I  )/NIXI 
EXTXR(  n-EXTXR(  D/NIXl 
505  CONTINUE 

DO  530  I -I .NT 
STirC-(I-l)*SNAP 
PRINT  560.STirC 

560  FORMAT!  I  HO.  •SNAP  TIME  -•.FIO.D 
J-MINO!NVX.It-INT!STIhE/VT)  ) 

PRINT  570. J 

570  FORMAT!  IHO.-VXLEYS  FIRED  -•.15) 

PRINT  540.NIX(I  )  .  (NSAMM!  I  .  IR)  .  IR-I  .NOEVIO 
5H0  FORMAT CO^. •NUMBER  TIMES  NO  SCREEN  •.I5.3!5X.*  IN  SECTOR  •.  15)) 
PRINT  600.  ( (SSL!  I  .  I  R)  .SSLS!  1  .  IRM  .  IR-1  .  NDEVIC) 

PRINT  601 .( (SSHd  .  IR)  .SSHS(  I .  IR) ) .  IR-l  .NOEVIO 
PRINT  602.((SLY!I.IR).SLYS(I.IRM.IR-I,ND£VIC) 

PRINT  603, ( (SRY! I . I R) ,SRYS( I , IR) ) . IR-I ,  NDEVIC) 

PRINT  604.EXTXL(1  )  .EXTXLSd) 

PRINT  605.EXTXRd  )  .EXTXRSd  ) 

600  FORMAT! •O^.*  SCREEN  LENGTH  •.2X,3!^  AVE  •.FB.  I,2X,^  SIGMA  •,F8. Id 

60)  FORMAT(»0^.»  INTER I  DR  HOLE  •,2X,3!»  AVE  •.F8.I.2X.^  SIGMA  • .F8. I ) ) 

602  FORMAT (•O*.^  LEFT  C OORO  •.2X,3!^  AVE  •,F8.I.2X,«  SIGMA  • .FS. 1) ) 

603  FORMAT(^0^.^  RIGHT  COORD  •.2X,3!^  AVE  •.F8.I.2X.*  SIGMA  • .F8. I ) ) 

60^  FORMAT! •O^.*  EXT  LEFT  C00R0^.2X,«  AVE  •.F8.l.2X.^  SIGMA  •.FOM 
605  FORMAT! •O^.^  EXT  ROHT  COOR0^.2X.^  AVE  •.r8.l.2X.*  SIGMA  •.F8.d 

530  CONTINUE 
PRINT  650 

650  FORMATdHO.  IBHOISTR  I8UTI0N  OF  SL  ) 

DO  680  l«) .NT 
STlME-d-n^SNAP 
PRINT  560. ST  I  rC 
J-MIN0(NVX.  I  +  INT(STIME/VT)  ) 

PRINT  570. J  ■ 

DO  678  IR-I .NDEVIC 
DO  675  L-I .4 
K2-25^L 
KI<K2-2H 

PRINT  660. rLEN(K) .K*K) .K2) 

660  FORMAT () HO . 25 1 5 ) 


126 


PAGC  16 


PRINT  670. IDSl  f I . IR ,K) »K«KI ,K2J 
670  rOPMATC IH0.2X,a5r5.  01 
PRINT  672 
872  FORMAT (I  HOI 
675  CONT I  NIX 
678  CONT  I  NIX 
680  CONT I  NUT 
RETURN 
END 
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SUBROUTlfC  STCLfQMUN, YIELD, SI0XS,<5IC.Y'i  ^lOZ^.ALPHAP.BEf AP.UPI 

COMMON  /V«:LD/  FSIGX,FSIGY,rsiOZ 
COMMON  /v«:ld/  dsigx.dsioy.dsioz 

COMMON  /yPCLD/  ASIOX.ASIGY, ASIOZ 
COMMON  /VPCIO/  ESIGX.ESIOY.ESIGZ 

COMMON  /v«:ld/  wndspo, factor 

WNOSPO-UP 

ES10X-G.9a94 

ES1GY»ALPHAP 

ESIOZ-BETAP 

FACTOR«QMUN»Y|ELD-. 156907 
ASIGX-DSIGX* ( (SIGXS/FSIGXI ( I .Q/ESIGXI  I 
AS1GY«D5IGY«C  CSIOYS/FSIGY»««(  I  .O/ESIOYI  ) 
ASIGZ»DSIGZ»I*SIGZS/FSIGZM»(I  ,0/ESI0ZI  I 

RETURN 

END 
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SUBROUriNE  FRMCLO(NT,OT,NVOL,NSl 


COmON  /WPSCR/  SIOX(501  ,SI0Y(50)  ,SIG7'*j01 
COMMON  /HPSCR/  FACIKSO) 


COMMON  /WPCLO/ 
COMMON  /WPCLO/ 
COMMON  /WPCLO/ 
COMMON  /WPCLO/ 
COMMON  /WPCLO/ 


FSIGX.rsiGY.FSIGZ 
0SIGX,0SI0Y,0SI02 
ASrOX.ASlOY.ASIOZ 
tSIGX.tSIGY.ESICZ 
WNOSPO, FACTOR 


PRINT  m 

!»♦  FORMATCO«.-NO.*,r5X.»Sr2E*l 
NS- 50 

00  199  l«[ .NS 
T-OTMl-ll 

IF  (T.GE.0.01  GO  TO  150 

srox(  n--i  .0 

GO  TO  199 
150  CONTINUE 

srox(  I  i-Fsrox*  (  (wnospo*t+asigx)  /osigxi  *  *esigx 

S1GY( I )»FSIOYM (WNOSPO'T+ASIOYl/OSIGY) • 'ESI  GY 
SrOZ(  1 1-FSr  GZ' ( (WNOSPO*  T+ASIOZl  /OSIGZ  I*  *ES  I GZ 
FACTTCI )»FACTOR/(SI 0X( 1 1 •SlOYf  n*SIOZ( I  1  1 
PRINT  16. 1 .SIGX(  II  .  SIGYC I) ,SIGZ(| I 
16  FORMAT!*  *. I4,6( IX, FIG.2, IXII 
199  CONTINUE 
RETURN 


END 
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SlJBROUTIf^  TIMFtMr.Ul 

COM^^ON  /REL/  XB160)  .YBteOJ  .2B<60I  .  TBf60  1  ,r«URST 
COMMON  /REL?/XCENTP (60) .YCENTP(601 

COMMON  /V^CR/  SlOX  (50)  .SIGYtSO)  ,S1G2(50) 

COMMON  /ABS/XBA(60) .YBA(60) ,2BA(60) .TBA (60) .NBRSTA 
COMMON  /ABS2/  XCENT (60 ) , YCENT (60 ) .2CENT (60 ) 

COMMON  /ABS3/  SIOXA (50) .SIGYA(50) 

COMMON  /TMSCLE/ INSNAP,  1NV0L,DT 
COMMON  /i#COtR/WVUI  .WVUB 
COMMON  /ACeiNO/  IAGE(60) 

DO  I  1-1 .NBURST 
AGE-(MT-I ) • INSNAP*DT-TB( I ) 

MUNGRP-AGE/DT+I 
1  AGE ( 1 ) -MUNGRP 
ir(MUNGRP.LE,0)  GO  TO  10 
XCENK I )-XBA( 1) ♦  U  'NVUl ‘AGE 
YCENT( 1 )-YBA( I )♦  U  •WVUa*AGE 
XCENTPC I )-XB( I )♦  U  ‘AGE 

S1GXA( n-AMAXl (ABS(SIOX(MUNGRP)»WVUI ) , AB5( SI  GY (MUNGRP) •WVUS ) ) 
SIGYA( n-AMAXl  (ABS(  S I GX ( MUNGRP )• WVUe ) , ABS(SIGY(MCWGRP) ‘NVUt  ) ) 
I  CONTINUE 
RETURN 
10  CONTINUE 
END 
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SUBROUTINE  DENSTYtXP. YP.2P.C0N) 

COMMON  /REL/  X0(6O1 ,Y0f6OI .ZB(60) .TBtaO J .NBURST 
COMMON  /R£:L2/XCENTP(60)  .YCENTPreO) 

COMMON  /WPSCR/  SIGX{50).SIGY(50),SIGZ(50J 
COMMON  /WPSCR/  TACTTISO) 

COMMON  /AGE  I  NO/  I  AGE (60) 

CON-0.0 

DO  399  MUN«=I  ,NBUR5T 

MUNGRP- 1  AGE  C  MUN ) 

IF(MUNGRP.LE.O)  GO  TO  400 

X=XP-XCENTP(MUN) 

Y=YP-YCENTP(MUN) 

Z=ZP-ZB(MUN) 

C0NN«(-.5*( (X/SIGX( MUNGRP) )**a 

*  ♦(Y/SIGY( MUNGRP) ) *'2 

•  ♦(Z/SIGZt  MUNGRP) )  ••2) ) 

IF(C0NN.LE.-25. )  GO  TO  399 
CON*CON+FACTT  ( MUNGRP )  *EXP(  CONN) 

399  CONTINUE 

*400  CONTINUE 
RETURN 
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SUBROUT ir^  LCONfPOX .POY.POZ.DX.DY,07,TOTl  NC>* 

COMMON  /REL/  XB(60) .YB(60I .ZB(60I . TB(60 1 »NBURST 
COMMON  /REL2/XCENTP(60I . YCENTP(eO> 

COMMON  /WPSCR/  S!GXC501.SIGYC50I,SIGZ(50) 

COMMON  /WPSCR/  FACT Tr 50) 

COMMON  /AGEINO/  IAGE<60) 

TOTLNC-0.0 

DO  399  N-I .NBURST 
I-IAGECN) 

IF( 1 .LE.O)  GO  TO  400 
ircSlGXM).LE.O.O)GO  TO  399 

A*(DX/SIGX( ! ) ) ••2*f  DY/S!GY( 1 ) ) • 'a* f OZ/S IGZ* I) ) ••2 
B'DX*(P0X“XCENTP(N) ) /SlGXf I) ••2+ 

•  DY- (POY-YCENTP(N) ) /S!GY( 1 ) ••2+ 

'  DZ'  (POZ-ZB(N)  )/S!  GZ(  1 )  "a 

S— B/A 

PMX-POX+S*DX 

PMY-POY+S'DY 

PMZ-POZ+S'DZ 

C-( (PMX-XCENTPfN)) /SlGXf ! ) ) ••2+ 

•  ( (PMY-YCENTP(N) )/SlGY( ! ) ) ••2+ 

•  ( (PMZ-ZB(N)  )/SlGZ  (  M  )  "a 
CC-C-.5-C) 

!F(CC.LE.“25. )  GO  TO  399 

FACLlN-FACTTf I ) -B.  50662B3*EXP(CC ) /SORT ( A) 
TOTLNC“FACL ! N+TOTLNC 

399  CONTINUE 

400  CONTINUE 
return 

END 
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PROGRAM  SMOKEHCdNPUT.OUTPUT.TAPEH- INPUT  1 


REAL  M_CX:X,M_OCY,Ti.OCZ 
REAL  NSIZX,NSIZY,NSIZ2 
REAL  NLOCXA.NLOCYA.M.OCZA 
REAL  N5IZXA«NStZYA,  NSIZZA 

COMMON  XB(  100)  .YB(  I  00)  .ZBMOOl  ,TB(  1001  ,  NBURST 
COMMON  FLOCXt 100) .ELOCYC 100) .ELOCZt 1001 

COMMON  rsizx(ioo) .rsiZY(ioo) .rsizzt loo) 

COMMON  NLOCXC100),NLOCY(I00),NLOCZ(100) 

COMMON  NSIZX(100).NS1ZY(I00),NSIZZ(100) 

COMMON  /ABS/XBAT 100 ) . YBA( 100) ,ZBA( 1 00) , TBA( 100) .NBRSTA 
COMMON  /ABS/ELOCXAC  100) .ELOCYAC lOO) .ELDCZAC  100) 

COMMON  /ABS/rsiZXA( 100) .PSIZYACIOO) .ESIZZAC 100) 

COMMON  /ABS/NLOCXAt  100) .NL0CYA( lOO) ,NLDCZA( 100) 

COMMON  /ABS/NSIZXA(  1 00) ,NSIZ YA( 100) .NSI ZZA( 100) 

COMMON  /PLACE/  SIGBR .SIGBO.SIOAR ,S1 GAD.  REL ,RDI SP 
COMMON  /PLACE/  XIDE  AL(6)  .  YIDEAL  (6)  ,ZIDE  AL  ( 6) 

COMMON  /VWDOIR/WVUI  .WVUB 
COMMON  /ACE I NO/  lAOE(IOO) 

COMMON  /ELDSIT/  CUT OEL .CUTQER.X INC 


COMMON  /HC5CR/  XCENTtSO, 150) 

COMMON  /HC5CR/  SIOX (20 . 150 ) .S1DY(20 , 1  50 ) . SI GZ (20 , 1 50) 
COMMON  /HC5CR/  FACT T 120, 1 50) , WnjFF (20 ) 


COMMON  /HCCLD/ 
COMMON  /HCCLD/ 
COMMON  /HCCLD/ 
COMMON  /HCCLD/ 
COMMON  /HCCLD/ 


FSIGX.FSIGY.FSIGZ 
DSIDX.DSIGY.DSIGZ 
ASIGX.ASIDY.ASIDZ 
ESIDX.ESIDY.ESIOZ 
WNDS PO .  FACTOR . T I NCR 


DITCNSION  TITLE(B) ,CATTN(3) ,THRES(3)  .CLTHRS(3) 

DATA  FSIGX/3.m/.FSlGY/3.Hl/.FS1GZ/l  .35/ 

DATA  DSIOX/IOO.O/ .DSIGY/ 100.0/ ,DSlGZ/20 .0/ 

DATA  ASlDX/0.0/ .ASl OY/0 . 0/ . ASlGZ/0 . 0/ 


READ(H.2)  TITLE 

2  FDRMAT(BAIO) 

PRINT  1. TITLE 

I  FORMAT! IHI .BAIO) 

READ(H,3l  U, ALPHA, BETA 

3  FORMAT(3F10.51 
PRINT  H.U. ALPHA. BETA 

4  F0RMAT(«0«.*  WIND  SPEED  '.FIO.S. 

*  •  ALPHA  •.FI0.5.*  BETA  '.FIO-S) 


READ(4.1I)  WVUI.WVU2 

11  FDRMAT(2F10.51 
PRINT  I2.WVUI,WVU2 

12  FORMAT! 1  HO. ‘WIND  VECTOR  DIRECTION  COSINES  X  -•.FIO.S.'Y  **.ri0.5) 

READtH.B)  YIELD, QMUN. BURNT 
6  FDRMAT(FI0.5.E10.5. F10.5) 

PRINT  10, YIELD. QMUN. BURNT 
10  F0RMAT(«0«.«  YIELD  •‘•.F10.5, 

•  •  QUAN  OF  MUN  ADJ  -•.E10.5.«  BURN  TIME  *«.riD.5) 

CALL  STCL(QMUN. YIELD. BURNT, ALPHA. BETA, U) 

READ(H.B)  NDEVIC 
B  FORMAT ( 15) 

READ(4. 15)  ( (CATTNt  IR) ,THRES( IR) ) ,  lR-1  .  NOE VIC) 
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15  FORMAT fDf 10.51 

DO  159  l-l,NOeviC 

CLTHRSl 1 )-ALOG( 1.0- THRESl 1 J) / ( -CATlNt I  1  i 
PRINT  20, 1 .CATTN* 1 )  . THRESt I  I .CLTHRSf 1 1 
159  CONTINUE 

?0  FORMAT ( *0' . 13. •  COEF  ATTN  ••.riO.5, 

1  •  THRESHOLD  -•,F10.2, 

?  •CL  THRES  -•.F10.21 

C  REAO  IN  PLACEMENT  INFORMATION 

READl^t.SO)  NRPV.NMPR 
50  F0RMAT(2I 10) 

PRINT  m, NRPV.NMPR 

m  FORMAT!  1  HO.  •NUMGER  OF  ROUNDS  PER  VaLEY  •.!'-/ 

1  IHO, -NUMBER  OF  MUNITIONS  PER  ROUND*.  IH/ 

2  IHO, •IDEAL  IMPACT  POINTS  •/ 

3  IHO.*  NO.  *. BX. •X* , I IX,*Y*, I  IX. •Z*) 

DO  2500  J- 1 . NRPV 

REAOm.55)  XIDEAL!  J)  .YIDEAL!  J)  .ZIDEALU  1 
55  FORMAT (i»Fl 0.5) 

PRINT  16.  J.  XI  DEAL!  J  1  .YIDEAL!  J)  .Z  IDEAL  (J) 

2500  CONTINUE 

16  FORMAT!*  *. 14.61 IX. F10.2. IXn 

REA0(H.2l)  NSAMP.NVa.NT.DT.OL 
21  F0RMAT(3I5.2F5.0) 

PRINT  23.NSAMP.NV0L  .NT.OT.DL 

23  FORMAT  M  HO.  *TOTAL  NUMBER  OF  SAMPLES  -*.  15,/ 

1  IX.*T0TAL  NUMBER.  OF  VOLLEYS  -*.15./ 

2  IX. -TOTAL  NUMBER  OF  TIMES  -*.  15./ 

3  IX. -TIME  INCREMENTS  -- .F5. 0 . *5EC  .  *  .  / 

4  IX. -LINE  INCREMENTS  -*.F5.0) 

REAO (4. 19)  VT.SNAP 
19  FORMAT(2F10.  I  ) 

PRINT  25. VT.SNAP 

25  FORMAT!  IHO, -TIME  BETICEN  vaLEYS  -*.,F10.l./ 

I  IX. -SNAP  TIME  INCREMENT  --.FlO.  1) 

READIH.B21  SlGBR.Sl  OBO.SIGAR.SIGAD.REL  .  ROI SP 
62  FORMAT (6ri 0.4) 

PRINT  24. SIGBR.SIGBO, S1GAR.SI0AD.REL.ro  ISP 

24  FORMAT! IHO. *BR  •* .F5. 0.5X. -0D  --.FB-O./ 

1  lX.-AR  -*.F5.0.5X.*A0  -*.F5.0./ 

2  IX.-REL  -*.F5.0./ 

3  lX.-RAOlAL  OISP  --,F5.0I 

NBRST A-NVOL - NRPV- NM  PR 

XOBS-0.0 

YOBS-0 . 0 

REA014.I5)  YTLHE.CUTOFL  .CUTOFR 
REA014.15)  XT.YT 
REA0(4,15)  XINC 
Z-l  .5 

PRINT  74, XINC 

74  FORMAT! IHO, *X- INCREMENT  - - . F5 . 1 ) 

PRINT  75.CUT0FL.CUT0FR 

75  FORMAT! IHO. *90  DEGREE  SECTOR  -•.F7.0.-  .  •.F7.0) 

PRINT  73.X0BS.Y0BS 

73  FORMAT! IHO. -OBSERVER  COORO  --.2F10.1) 

PRINT  71 .XT.YT 

71  FORMAT! IHO.-AIMPOINT  COORD  *-.2F10.1) 

PRINT  72.YTLIfC 
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It*  FORMAT  ( 1  HO,  •DISTANCE  TO  L  I  NE  '•,F!0  M 
PRINT  35, Z 

35  FORMAT  (*0*, 'HEIGHT  ABOVE  CCNTROIO  51 

CALI  CLEAR  (DL) 

DO  1000  KI«I .NSAMP 
PRINT  7.KI 

7  FORMAT! IHI , 'DATA  FOR  SAMPLE  NUMBER  ',151 

CALL  MPLACE(NVOL,NRPV,NMPR,VT,XT,YTI 

CALL  CONBST 

PRINT 

FORMAT!  •O'.'  MUNITION  BURSTS  AT  •) 

PRINT  M3 

M3  FORMAT! 'O'  ,  I  I  X. 'LOG  AT  ION' ,  I  I  X.5X.7X  .  •  T  I  ME'  ) 

1=0 

DO  199  IV- I ,NVOL 
DO  199  IRO-I ,NRPV 
DO  199  1 5- 1 .NMPR 
I-I^^I 

PRINT  MB. I , IV, IRO,  I  S.XBA! I  1  .YBA! I )  ,ZBA(  I ) , TBA( I) 

199  CONTINUE 

MB  FORMAT!'  • .M ( I  X. 1 3 ,  I X) , 3( I  X. FB.B. I X) .6X  ,  FB-B) 

DO  399  IT- I ,NT 
T-FLOAT!  I  D'SNAP 

PRINT  ae.T  . 

?e  FORMAT! 'O', 'TIME  AFTER  EMISSION  '.FIO.O) 

CALL  TIMEITI 
CALL  CON 

C  PRINT  32 

32  FORMAT!'0'.'  MUNI T I  ON  CHARACTERI ST  ICS  ') 

C  PRINT  33 

33  FORMAT! 'O' . I  I X. 'LOC AT  ION' . I  I X .5X. IMX, *S IZE' 1 
DO  350  IMUN-I .NBRSTA 

C  PRINT  3B, IMUN.NLOCXA! IMUN) .NLOCYA! IMUN) .NLOCZA! IMUN) , 

C  •  NSIZXA!  IftJN)  .NSIZYAI  IMUNJ  .NSIZZA!  IMUN) 

C  PRINT  39. IMUN.FLOCXA! IMUN) .FLOCYA! IMUN) .FLOCZA! IMUN) , 

C  •  FSIZXA!  IMUN)  .FSIZYA!  IMUN)  .FSIZZA!  IMUN) 

350  CONTI NIC 

3B  FORMAT! 'O' . IX.  12, IX . 'hCAR* . IX,3! I X,FB.2 , IX) .5X. 3! IX.FB.2. I  X )  ) 
39  FORMAT!'0'.IX.I2.IX.'  FAR' . I X . 3 ! I X . FB . 2 . 1 X ) . 5X . 3 ! I X , FB . 2 . I X ) 1 

CALL  SIZE!XLOW.XHIGH.YLOW. YHIOH) 

CALL  MXMIN!YTLIh€:.XMIN.XMAX) 

PRINT  5.XMIN.XMAX 

5  FORMAT!IH  .'XMIN  -•  .FI  0 .2.2X . 'XMAX  -'.FIO.B) 

YSLIhC-500.0 

C  CALL  MATCON!Z.XLOW. XHI GH , YLOW, YHI GH) 

C  CALL  MATCL!Z.XLOW.XHIGH.YOBS.YSLINE) 

CALL  CALC  !XOBS.  YOBS.  Z.YTLIfC.XMIN.  XMAX.  NOEVIC.CLTHRS) 

CALL  ENDPTS!XMIN.XMAX. IT) 

DO  299  IR«I .NDEVIC 
CALL  VEVAL  ! IT. IR) 
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299  CONTINUE 
399  CONTINUE 
1000  CONTINU 

CALL  CALPRT  f  N9AMP,  NT  .SNAP.NVa  ,  VF  .NOf  V  IC » 

STOP 

ENO 
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SUBROUTINE  CLEARfDL) 

COMMON  /STATS/  SSL  I  50 . 3)  . SSLS( 50 , 31  , SI  Y  ( 50 . 31  , SL  YS( 50 . 3)  , 

1  SRYC50. 3) .SRYSfSO, 3l .SSHI 50 , 3) , SSHSI50 . 3» .NSAMM(50. 3) .LENf  1001 . 
a  DSLf50.3, 1001 

COMMON  /STATa/  EXTX L C 50 1 . EXTXL S ( 50 » , EXT XRI 50 1 . EX FXRSf 50 1 . N I  X C 50 ) 

DO  800  T-l , 100 
800  LENI I ) -DL* I 
DO  laO  l-l .50 
NIXCn-0 
EXTXL(  n-0.0 
EXTXLSI I )-0.0 
EXTXRI I »-0.0 
EXTXRSI  n-0.0 
DO  no  IR-I  .3 
NSAMMn,lR»-0 
S5L( I . |R)-0.0 
SSLS( I . IR)-0.0 
SLY( I . IR)-0.0 
SLYSI I . IR)-0.0 
SRYd  .|R)-0.0 
SRYS( I . lR)-0.0 
DO  105  K-l .100 
105  DSL(  I  nR.K)-0.0 
SSHI I .IR)-0.0 
I  10  SSHS( I . IRl-0.0 
laO  CONTINUE 
RETURN 

END 
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28 

30 

32 


SU0ROUT|^C  MPLACEfNVOL.NRPV.NMPR.VI.vt  .  YTI 

COHHON  /ABS/  XB(  1 001  ,  YB(  I  001  ,ZB(I  OU  <  .  T  B(|  001  .^eURST 
COMMON  /PLACE/  5IGBR,5I  OBO.SIGAR.SIGAO,  REL  .ROISP 
COMMON  /PLACE/  XIDE  AL  ( 6) , YIOEAL (61 .ZIOE  AL (61 
CALL  N0RAN(R.5IGAR,E.5IGA0) 

XCA-XT+E 
YCA-YT^R 
DO  32  I -I ,NVOL 
0RSTTM-  ( I-I 1*VT 
DO  30  J-I ,NRPV 

IFCRANr(OUM) .OT.REL  1  GO  TO  30 
CALL  NORANIR.SIGBR.E.SIGBDI 
XCRD-XCA^XIOEALI  JHE 
YCRO»YC  A^  Y I  DEAL  ( J1  ♦  R 
DO  28  K-l ,NMPR 

MUNNO»(  I-II»NRPV»NMPR  ♦(J>n*NMPR  ♦K 
CALL  MUNOSP(ROISP.XMUN,Y^OJ) 

XB(HUNNO)  •XCRDi^XMUN 

YB( MUNNO) -YCRD^ YMUN 

ZB(MUNN01-0.0 

TBCMUNNOl -0RSTTM 

CONTINUE 

CONTINUE 

CONTINUE 

RETURN 

END 
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SueROUT I  CON0ST 

COMMON  XB(  !00)  ,  YBf  I  OOJ  .ZBMOOl  ,  TBI  lOni  .  NOURST 

COMMON  /ABS/XBAf  100  I.,YBA<  lOOt  .ZBA(  100)  ,  TBAI  100)  .NBPSTA 
COMMON  /wNooiR/wvui  , wvua 

N8URST«N0RSTA 

DO  100  I-I .NBRSTA 

XB< I )*WVUI 'XBAC I )+WVU2*YBA( I ) 

YB< I )«-WVUa*XBA(  I ) ♦WVUI 'YRAf I ) 

ZB<n*7BAU) 
rB» I )*TBA< I  1 
100  CONTINUE 
RETLJRN 
CNU 
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SUBROUTINE  CON 

REAL  M.OCX,NLOCY,NLOCZ 
REAL  NSIZX,NSIZY,NS17Z 
REAL  NLOCXA.NLOCYA, NLOCZA 
REAL  NSIZXA,NSIZYA,NSIZZA 

C0^f10N  XB(  100)  ,YB(  I  00)  ,ZB(  I  OOl  ,TB(  100) ,  N0UR5T 
co^f10N  rLocx(ioo).rLocY(ioO),rLocz(ioo) 
co^f10N  rsizx(  1001  .rsizYdoo)  .rsizzc  loo) 

COrtlON  M-OCX(|QOI,NLOCY(IQO).NLOCZ(IQQ) 

COrtlON  N5IZXU00).NSIZY(I00).N5IZZ(I00) 

COMMON  /ABS/XBA( 100 ) ,YBA( 100) *ZBA( 100) ,  T8A( 100) .N0RSTA 
COMMON  /ABS/ELOCXAC  100) .ELOCYAC 100) ,ELQCZA( 100) 

COftlON  /ABS/r5IZXA(  1 00) ’.ESIZYAC  I  00 1  ,ESI  ZZA(  lOO) 

COftlON  /ABS/M_OCXA(  100)  ,M_OCYA(  100)  ,M-QCZA(  lOQ) 

COMMON  /ABS/N5IZXA( 100) ,N5IZYA( 100) ,NSI ZZA( 100) 

COMMON  /V#JOOIR/WVUl  ,WVUB 

DO  100  I-I,NBUR5T 

NLOCXAC  I  )-WVUl  'M-OCXt  I  )-WVU2*M.0CY(  I  ) 
hLOCYAd  )-WVUe*NLOCX(  I  )  ♦WVUI  •^^.OCYd  ) 

NLOCZAd  )-NLOCZd  ) 

NSIZXAC I )-NSIZX( I ) 

NSIZYAC  |)-NSIZYd) 

NSIZZACI  >-NStZZd  ) 

ELOCXAt  I  )-WVUI*ELOCXd  )-WVUa*rLOCY(  I  ) 

rLOCYAd)-Nvua«rLocxd)+wvui  •elocy(  i  )  ■ 

FLOCZACD-FLOCZd  I 
FSIZXAd)-FSIZXd  ) 

FSIZYACD-FSIZYl  I  )■ 

FSIZZAl I )-FSIZZ( I ) 

100  CONTINUE 
RETURN 
END 
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SUBROUTirC  MUNOSPIR  ,XMUN,YMUN) 

RS-SQRT  I  RAhr  ( DUM  M  •  R 
THETA-RANTIDUMJ -B.PBaiB 
XMUN*RS*COS( ThCTAI 
YMUN*RS-SINCTh€TA) 

ReruRN 

END 
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SLIBROUriNE:  NORAN(R.  SR.D.SO) 
I  X»RANF(DUM1 

irix.LE.o.OJ  00  ro  i 

A*SQRTl -e.O'ALOGlXl  \ 

B*6 . 2B3 1 85307 1 0*RANF I DUM I 
R*A*SR*SINrB» 

D«A*S0*C0S(BI 

RTTURN 

FNO  ■ 
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SUBROUTIhC  SIZEeXLOW.XHIGH.YLOW,YMIGH» 

REAL  M.0CX,M.0Cy,M.0CZ 
real  N5IZX,N5IZY,N5IZZ 

COMON  /ABS/XBf  100)  ,  YB(  I  00 ) ,  ZB(  1 00 ) ,  TB(  100)  ,N0URST 
COMON  ./ABS/FL0CX(  I  00)  ,rL0CY(  100)  ,FL0CZ  ( 100) 

COMMON  /ABS/rsiZX( I  00) .FSIZYf 1 00 ) ,FSIZZ  C 100) 

COMMON  /ABS/NLOCX(  I  00) ,M.OCYt 100) ,NLOCZ ( 100) 

COMMON  /AB5/N5IZX( I  00) ,NSIZY( 100) ,N5IZZ  MOO) 

XLOW-I.OE+6 

XHIGH--I .OE+6 

YLOW«l .OE+6 

YHIOH— l,0E+B 

DO  10  I-I,NBURST 

TEST-CNLOCX( I)-H.«NSIZX( I ) ) 

IF(TEST.LT.XLOW)  XLOW-TEST 
rF(TEST.GT.XHIOH)  XHIGH-TEST 
TEST-(NLOCX(  I  )'m.«NSIZXM  )  ) 

IF(TEST.LT.XLOW)  XLOW-TEST 
ir(TEST.GT.XHIOH)  XHIGH-TEST 
TEST-(FLOCXM)-H.TSIZXM) ) 

IFCTEST.lt. XLOW)  XLOW-TEST 
IFCTEST.OT.XHIGH)  XHIGH-TEST 
TEST-CFLOCXC I )^H.«rSIZX( I ) ) 

IF(TEST.LT.XLOW)  xlow-tcst 
IFCTEST.OT.XHIGH)  XHIGH-TEST 
TEST-CNLOCYM  I-^.^NSIZYC  I  )  I 
IFtTEST.LT.YLOW)  YLOW-TEST 
IFCTEST.GT. YHIOH)  YHIGH-TEST 
TEST-CNLOCYC I ) +H . -NSIZYC I ) ) 

IFCTEST.LT. YLOW)  YLOW-TEST 
IFCTEST.GT. YHIOH)  YHIGH-TEST 
TEST-CFLOCYC  I)-4.*FSIZYM)) 

IF CTEST.lt. YLOW)  YLOW-TEST 
IFCTEST.GT. Y)^IOH)  YHIGH-TEST 
TEST-CFLOCYC I )*4.«FSIZYC I ) ) 

IF CTEST.lt. YLOW)  YLOW-TEST 
IFCTEST.GT. YHIGH)  YHIGH-TEST 
10  CONTINUE 
RETURN 
END 
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SUBROUTINE  MXMIN( VTL 1 1C  . XMIN.XMAXI 

REAL  NLOCX.NLOCY.NLOCZ 
REAL  NSIZX.NSIZY.NS IZZ 

COMMON  /ABS/XBI 100) ,YB( 100) ,ZB< 100) ,TB( 1001 .NBURST 
COMMON  /ABS/rLOCX(  I  001 .FLOCYt 1001 .FLOCZt  1001 
COMMON  /ABS/FSIZXl  I  00) ,F5IZY( 1 00 ) .FSIZZ (  1001 
COMMON  /ABS/NLOCX(  I  00) ,NLOCY( 100) ,NLOCZ  C  1001 
C0^f10N  /ABS/NS1ZX(  I  00)  .NSIZYC  1001  ,NSIZZ  (  1001 

RATMX"-I .OE+IO 
RATMN-I .OE+IO 
DO  10  1-1 , NBURST 
1F(NL0CY( I) .EQ.0.0)  GO  TO  5 
RATL-(NLOCXC I )-^.*N5IZXri 1 )/NLOCYt 1 1 
IF(RATL.LT.RATht<)  RATMN-RATL 
IF(RATL.OT.RATMX)  RATMX-RATL 
RATN-tNLOCX( I )+H.*N5IZX( I ) )/NLOCYI I  1 
IF(RATN.LT.RATW)  RATMN-RATN 
IF(RATN.OT.RATMX)  RATMX-RATN 
5  CONTINUE 

IFtFLOCYl I ) .EQ.0.0)  GOTO  10 
RATL-(FLOCX(  I ) -»♦.  •FSI  ZX(  I )  )/FLOCYC  I  1 
IF(RATL.LT.RATMN)  RATMN-RATL 
IF(RATL.OT.RATMX)  RATMX-RATL 
RATN-(FLOCX( I  1 . •F 51 ZX ( I ) 1 /FLOCYl 1 1 
IF(RATN.LT.RATMN)  RATMN-RATN 
IF(RATN.OT.RATMX)  RATMX-RATN 
10  CONTINUE 

XMIN-YTLIIC^RATMN 

XMAX-YTLIIC»RATMX 

RETURN 

END 
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SLOROUTINE:  CALCCXOBS,  VOBS,Z,YTLirC  .  .XMAX.NOFVIC.CLTHRS) 

COMMON  /FLDSIT/  CUTOTL .CUTOTR^XINC 
COMMON  /LINES/  V( 1001 ,31 , IXLEFT, IXRGHI 
COMMON  /WNOOIR/WVUI  ,WVU2 
DlfCNSION  CLTNRS(3i 

XL-AMAXl (CUTOFL.XMI  N) 

XR-AMINI (CUTOFR.XMAXI 

IXLEFT-INT( (XL-CUTOFL)/XINCI  ♦  I 

|XR0HT-|NTC(XR-CUTDFL)/XINC)  ♦  I 

0031  IP-IXLEFT, IXRGHT 

X-XINC  •  FLOAT! IP-1)  +CUTOFL 

OX-X-XO0S 

OY-YTLIlE-YOeS 

DZ-0.0 

VMAG-5QRT(0X»»2+0Y* •2+0Z»»2) 

OX-OX/ VMAO 
OY-OY/VMAG 
OZ -07/ VMAO 
OXP-WVUI •0X^WVU2»0Y 
0YP--WVU2*0XfWVUI •OY 

CALL  LCONlXOBS.YoeS ,Z.OXP.OYP.OZ.TOTLNC  1 
00  29  IR-I .NOeVlC 
V( IP, IR)-0.0 

IFCTOTLNC.OE.CLTKRS!  1R1 1  V( IP, 1R1-1 .0 
29  continue 
31  CONTINUE 

00  35  IR-I .NOEVIC 

C  PRINT  3H. (V( 1 . IR) . I -IXLEFT. IXRGHT) 

35  CONTINUE 

3^  FORMAT ( 1  HO . ( 25FH . 1 . / 1 1 
RETURN 
END 
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SUBROUTINE  VEVALCII.IR) 

COMMON  /ELDSIT/  cutotl.cutotr.xinc 

COMMON  /LlfCS/  V{  1001 ,31  ,  IXLEFT,  IXRGHT 

COMMON  /STATS/  SSL  I  50 , 3 1 . SSLS( 50 . 3 1 , SLY ( 50 . 3) . SL Y5C 50 . 3 1 , 

1  SRY(50,31 ,SRYS(50. 3) ,SSH(50,31 .SSMStSO ,3) ,NSAMMt50.3) .LEN( 1001 , 
a  OSL(50.3, 1001 

INTEGER  RY.SH 

LY«0 

RY«0 

SH-0 

SL  *=0.0 

LEFT  Y'COORO  «  LY 

IF( IXRGHT-IXLEFT.LE .01  GO  TO  95 
00  10  l-IXLEFT. IXRGHT 
IFIVC 1 . !R1 .EO.0.01  GO  TO  lO 
LY'XINC*FL0AT(  !“!  1  ♦CUTOTL 
K«I 

GO  TO  15 
10  CONTINUE 
GO  TO  HO 

RIGHT  Y-COORO  -  RY 

15  00  ao  J«K, IXRGHT 

IF(V( J. IRl .EQ.0.01  GO  TO  aO 
RY-XINC-FLOATC J-I  1  ♦  CUTOTL 
M-J 

aO  CONTINUE 

TOTAL  SCREEN  LENGTH  «  SL 
SL-RY-LY 

IFISL.EQ.O. 1  GO  TO  MO 
SSL  I  1 1 , IRl-SSLl  1 1 , 1  R) ♦SL 
SSLSni.IRl-SSLSlII  ,IR1^SL*SL 
SLY( I1,IR1-SLYCII,IR)+LY 
SLYS(Il,!R)-SLYS(I! .IRl+LY-LY 
SRYl 1 1 . IR1-SRY( I ! ,  I  Rl+RY 
SRYS( I  1 . IRl-SRYSt I  1  . IR)+RY*RY 


INTERIOR  HOLE  LENGTH 


SH 


DO  30 

IFCVC I , IRl .EQ. I .01  GO  TO  30 
SH-SH+XINC 
30  CONTINUE 

SSH( 1 1 . IRl -SSHl II .  I  Rl+SH 
SSHSI I  1 , IRl-SSHSl I  I  . 1 Rl+SH* SH 
GO  TO  50 

NSAMM  TO  BE  USEO  TO  MODIFY  NSAMP  SUCH  THAT 
AVE  LY,RY,ANO  SH  WILL  BE  CALC.  GIVEN 
NON-2ERO  SCREEN  LENGTH. 

MO  NSAMMI 1  I . IRl-NSAMMl  I  I , lRl  +  1 

DISTRIBUTION  OF  SCREEN  LENGTH 

50  00  60  I  - 1  . 1 00 

IFCSL.GE.LENt Ill  GO  TO  60 
OSL( I  I , IR.I 1-DSLl I  I  . IR. I l  +  I .0 
GO  TO  70 
60  CONTINUE 
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70  CONTINUE 
RETURN 
95  CONTINUE 

N5AW1f  I  I  .  IR)  «N5AMM(  I  I  ,  IR)  ♦  | 
Rf  TURN 
f  NO 
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SUBROUTIfC  ENOPTSCXMIN.XMAX, IT) 

COMMON  /STAT?/  EXTXL(50I .EXTXLS(501 .FXl XRJ50I ,EXTXRS(50) .Nl X(50J 
IFCXMAX-XMIN.LE. I .0 )  GO  TO  05 
EXTXL( IT)-FXTXL( IT)  ♦  XMIN 
EXTXLSMT)*EXTXLSt  I  T  )  ♦  XMIN*XMIN 
EXTXR(IT)«EXTXR(IT»  ♦  XMAX 
EXTXRSC IT)-EXTXRS( I T)  ♦  XMAX'XMAX 
RETURN 
05  CONTINUE 

NIXl IT)-NIX( I T»  *  I 

RETURN 

END 
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SUBROUTINE  MATCONt  Z  ,XL0W,XH10H,  YLOW,  YHI  GH) 

COMMON  /WNOOIR/WVUI  ,WVU2 
DIMENSION  XCOROmO)  .YCOROmO) 

DlfCNSION  CONMATr40,401 

DO  1 100- 1-1 *40 

XCOROI 1 ) -FLOAT  11 ) • ( I . 0/40 . 0 ) • ( XHIOH-XL  OW) ♦XLOW 
1100  CONTINUE 

DO  1500  l-I  .40 

YCOROI  n-FLDATl  I-l  1  •(  I  .  0/40 . 0)  *  t  YHIGH- YLOW)  ♦  YLOW- 
1500  CONTINUE 

PRINT  BB 

BB  FORMAT  (*>  X  '.SX.BOX,*  Y 

PRINT  B6. ( YCORD( IP) . IP-I ,B0) 

B6  FORMAT(«0'.7X.B0(F6. n  ) 

DO  B400  I-I .40 
DO  B400  J-I .40 

X-WVUl  'XCOROC 1 )  ♦WVUB«YCORO(  J) 

Y--WVUB»XCORO(  1  )+WVUI  •YCOROt  J) 

CALL  OENSTYCX.Y.Z.CON) 

CONMATCI ,J)-CON 
B400  CONTINUE 

DO  B500  1-1.40 

PRINT  B8.XC0R0( 1 ) , ( CONMATC I .IP) , IP-1 , BO  ) 
IBLANK-5*! 1/5) -I 
IF  IIBLANK.EQ.OIPRINT  3B 
B500  CONTINUE 
PRINT  BB 

PRINT  B6.(YCORO(IPJ  . IP-Bl .40) 

DO  3000  1-1,40 

PRINT  B8,XC0R0( I ) .  (  CONMAT( I . IP) , IP-BI .4  0) 
lBLANK-5MI/5)-I 
IF  (IBLANK.EQ.O) PRINT  3B 
3000  CONTINUE 

B8  FORMAT (•  • , IX.FB. I . B0(F6.0) ) 

3B  FORMAT (•  •) 

RETURN 

END 
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SUBROUTIIC  MATCL<Z,  XLOW.XHIOH, YLOW. 

COMMON  /WNDOIR/WVUI  .WVU? 

DIMENSION  XI  (‘♦01  .X^  mOI 
DIMErJSION  CONLINC^O  .HO) 

YI*YLON  • 

Y2*YHI(5H 

DZ*0.0 

DO  I  100  I-I  ,H0 

XI ( I ) -FLOAT ( I -I )•( I .0/H0.01-(XH1GH-XLOW)^XLOW 
I  100  CONTINUE 

DO  1500  I-I ,H0 

X2( I  1 -FLOAT ( I  - n • (  I  0/H0,0) • (XHIGH-XLOW 1 ♦XL ON 
1500  CONTINUE 

PRINT  I5,Y1.Y2 

15  FORMAT(«0«.«  YI  -.Fa.S,*  Y2  VEa.S) 

PRINT  BB 

BB  FORMAT  I •  XI  •.5X,20X,-  xa 

DO  5500  I-I .HO 
DO  2H00  J-I .HO 

DIST-SQRT(  (Xai Jl-X  U I  1 )-«a*(Y2-Yl )**ai 

DX-(X2(J)-XI ( I n /DI ST 

DY-(Y2-YI l/DIST 

X-NVUl *XI ( I ) ♦WVUe«YLOW 

Y--WVie«XI  (  I  1+NVUM  YLOW 

oxp-wvui •ox+wvue«DY 

DYP--NVIC«DX+NVUI "DY 

CALL  LCON(X,Y.Z,DXP,DYP.DZ,CON) 

COr<.IN(  I  ,  J)-CON 
anoo  CONTINUE 
2500  CONTINUE 

PRINT  26, (X2(I) ,1-1  .201 
DO  29  K-l .HO 

PRINT  20. XI (K) . (CONLINIK, 11,1-1,201 
29  CONTINUE 

PRINT  I5,YI,Y2 
PRINT  22 

PRINT  26. (X2l 11,1-21  ,H0) 

DO  07  K-I.HO 

PRINT  aa.XI (K) , (CONLINIK. 11,1-21 ,H01 

27  CONTINUE 

26  format t«0«.7X. 20 IF6. 1 11 

28  FORMAT (•  • ,F6. I .20( F6.0) 1 
RETURN 

END 
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SUeROUTIhC  CALP«TINSA#^,NT,SNAP,Mva.,vr  .NOCVICI 

■  C  CALCULATE  AlO  PWINT  STATISTICS 

COM-CN  /STATS/  SSLt  30,31  ,SSLSI50,3J  .SLY«50,3l  ,SLYSf30.3» . 

I  SRYI  50.31 .5RYSI30,  31  ,SSHf30.3I  ,SS»^(30 .31  ,NSA^rH30,3I  .LENI  1001  . 
a  0SL(30.3.I00I 

COmON  /STAT2/  EXTXLI30I .EXTXLSI30) ,EXTXRI30J .EXTXRSI30) ,Nl XI30J 
PRINT  SIO.NSAMR 

310  rORMAT<lHI,*TOTA*.  NUTCER  OT  SA»«^ES  IS* 

00  503  1-1 .NT 
DO  320  IR«I .NOeviC 
NSAf^l  -N!5A»1^  >ISA>^T  I  .  I  Rl 
IFINSA>^I  .L£.  )  IGO  TO  315 

SLYSI I  .  IRI-SORTHSLYSI  I  .  IRI-SLYI  I .  IR»  ••2,/NSA#^I  )/INSAMPl-!  1  ) 

SRYSI  I  .  IRI-SORTC ISRYSI I  .  IR)-SRYI  I .  IR)  •  •  2. /NSAr»»n  /  « NSA^f»l  - 1  1  ) 

S5HSI  I .  |R)-SORTUSSMSI  I  .  IR)-SSHI  I  .  IRl  ••  2. /NSAl-*»l  »/ f  NSAI-PI  - 1  l  ) 

SLYIl . IR)-SLYC I . IR) /NSAMPt  , 

SRYI  I .  IR)-SRY(  I .  IRI  /NSAJhPI 

SSWI  I  .  IR)-SSHI  I .  IR)  /NSAM»I  '• 

313  SSLS!  I  .  IRI-SORTC  ISSLSM  .  IRl-SSLI  I  .  IRI  ••2, /NSA«>)  /  (NSAf^-l  J  I 
320  SSLII.IRI-SSLCl.lRI/NSAf^ 

NIXI-N5Af«*^IX(  1 1 
iriNIXI  .L£.  1 1  00  TO  303 

EXTXLSC  M -SORT  U EXT  XLSI  M-€XTXL(  I  I  ••2- /  NI X!  I  /  INI  XI -I  J  I 
CXTXRSC I I-S0HTUEXTXR51 1 1 -CXTXRI  n  ••2.  /  Nixn  /  INIXI -m 
EXTXL 1 1 1-CXTXL 1 1 1 /N  IXI 
EXTXRI I I-EXTXRI  1 1 /N  I XI 
503  CONTINUE 

DO  330  1-1 .NT 
STIME-I|J*SNAP 
PRINT  SSO.STIfC 

560  rOFWATUHO.*SNAP  TlfC  -•.FIO.II 
J-HINOINVOL,I*INTtST|hC/VT)  » 

PRINT  570. J 

370  rORHATU HO. •volleys  FIREO  --.ISI 

PRINT  5^0.NIX(  U .  INSAWMI I . IRI . IR-I .NOEVIC) 

5*0  FORMAT  U  MO. -NmaER  TIhCS  NO  SCREEN  •  .  1 3 . 3*  2X.  •  IN  SECTOR  *.13)  > 
PRINT  600. 1 ISSLI  I .  I  R) .SSLSI I . IR* » . IR-I  .  NOEVIC) 

PRINT  601  .  (ISSHI 1  .  IR)  .S5HS(  1 .  IR)  1 .  IR-I  .NOEVIC) 

PRINT  602. 1 1  SLY  1 1. IRI .SLYSI I, IR)) .IR-I .NOEVIC) 

PRINT  603, II SRYI 1,1  R).SRTSI I, IRI), IR-I,  NOEVIC) 

PRINT  6CV.EXTXLI1),CXTXLSI1) 

PRINT  603, EXTXRI  I )  .  EXTXRSI 1 1 

600  FORMATI*0*.*  SCREEN  LENGTH  -.2X.3I*  AVE  *,F8.I.2X.*  SIGMA  <,ra,I)) 

601  FORMAT  I -O*.*  INTERIOR  HOLE  '.ax.  31*  AVe  •.F9. 1  .2X.*  SIGMA  •  .F9.  I  )  ) 

602  FORMATI-0*.*  left  COOflO  •.2X.3I*  AVE  •.F9.  I.2X.*  SIGMA  •.ra.M) 

603  FORMATI-0*.*  right  COOflO  •.2X.3I*  AVE  •.F9.I.2X.*  SIGMA  •.r9.I)) 

60N  FORMATI-0*,*  EXT  LEFT  COORD*. 2X.*  AVE  •.F9. 1.2X.*  SIGMA  -.Fa;!) 

603  FORMAT! -o*.*  EXT  ROMT  COOflO*. ax.*  AVE  -  .Fa. I.2X.*  SIGMA  '.rs,  M 
530  CONTINUE 

PRINT  630 

530  FORMATI  IHO.  l8)CISTfl  IBUTION  or  SL) 

DO  690  I -I .NT 
5TI^E-I I)*SNAP 
PRINT  360.5riT€ 

J-MINOINVOL.  I-INTIST|T€/VT)  I 

PRINT  570. J 

DO  679  IR-I. NOEVIC 

DO  673  L-I.V- 

K2-23*L 

KI-K2-eN 

PRINT  660.  ILENIK)  ,K-KT  ,K2) 

660  FORMAT  I  IH0,a3l3) 

PRINT  670.l0SLII.IR.K).K-KI,K2). 

670  FORMAT!  IH0.2X.25F3.  0) 
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6  rp  T  ORMA  T  t  1  HO  > 
P.7*)  CONTI  NIC 
6  70  (ONTINU 
680  CONTI  NIC 
FC  TURN 
TNO 
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SUBRCXJTirC  STCL CQMUN, YIELD, BTIME,AIPHAP,BCTAP,UPI 

cort«N  /HCCLD/  rsiox,rsiov.rsi02 
COmON  /HCCLD/  DSI0X,DS10Y,DSI0Z 
COmON  /HCCLD/  ASIOX,A5IOY,A5IOZ 
COmON  /HCCLD/  E5IOX,E510Y,E5IOZ 
COmON  /HCCLD/  WN05P0, FACTOR, TINCR 

^f>UFTL■I^O 
PUTT-FLOAT  (r«VTL » 

VWOSPO-UP 

ESIOX-ALPHAP 

ESIOY-ALPHAP 

ESIGZ-BETAP 

FACTOR-QMUN'YIELD*.  IP69tf7/PUFT 

TINCR-BTIME/PUFT 

RETURN 

END 
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SUeROUTI^€  TirC(T) 


REAL  ^^_OCX,M.OCY,^LOCZ 
REAL  NSIZX.NSIZY.NSIZZ 


COMMON  XB(  lOOJ  ,YB(  I  00 ) . ZB( I  00 ) , TB(  1 00 )  ,  NBURST 
COMMON  FLOCXC lOOJ ,FLOCY( lOOJ ,rLOCZ( too) 

COMMON  FSIZXI 100) .FSlZYt 100) ,FSIZZ( 100) 

COMMON  M.OCX(IOO),NLOCY(IOO),NLOCZ(tOO) 

COMMON  NSIZX(IOO)  .NSIZYdOOI  ,NSIZZ(I00) 

COMMON  /AOCtNO/  lAOE(lOO) 


COMMON  /HCSCR/  XCENKPO,  150) 

COMMON  /HCSCR/  SIGXtPO, 1501, SrOY(20, 150), SIGZ(P0, 150) 
COMMON  /HCSCR/  FACT T (PO ,I 50 ) ,NPUFF (20 ) 


COMMON  /HCCLO/ 
COMMON  /HCCLO/ 
COMMON  /HCCLO/ 
COMMON  /HCCLO/ 
COMMON  /HCCLO/ 


FSIOX,FSIOY,FSIOZ 
DSIOX,DSIGY,DSIOZ 
ASIOX.ASIOY.ASIGZ 
ESIOX.ESIGY.ESIOZ 
WNOS  PO ,  F  AC  TOR ,  T I  NCR 


PAGE-*1000. 

HJNGRP-0 


DO  299  MUN- I , NBURST 

AGE-T-TB(MUN) 

IFtAGE.I^.PAGE)  GO  TO  150 
IA0E(MUN)-MUNGRP  • 

GO  TO  205 

150  CONTINUE 

MUNGRP-MUNGRP+ I 

IAOE(MUNl*MUNORP 

PAOE-AGE 

DO  199  [-1 . 120 
TW-AOE-TINCR«FLOAT( l-l 1 
IFtTW.LE. 0.0)00  TO  200 
XCENT  ( MUNGRP ,  1 1 -WNOSPO* TW 

SIOX(MUNGRP.I)-FSIOX«t(WNOSPO«TW+AStOX) /DSIOXl ••ESIGX 
SIOY(MUNORP.I  )-FSIGY«((l#CSPO«TW+AS!OY)  /DSIOY)  ••EStCY 
SIOZ(MUNO«P,I)-FSIOZM(W«)SI>0*TW+AS!OZ)  /DSIOZI  ••ESIGZ 
FACTT(MUNORP,l )- 

1  FACTOR  /  (SI GX(MUNGRP.  I  )  -SIO  YtMUNGRP.  I  )  *51 GZ  t^«JNGRP,  I  )  ) 

199  CONTINLE 

NPUFF(MUNGRP)-120 
GO  TO  210 

200  CONTINLE 
NPUFF(MUNGRP)-I-I 


205  CONTINLE 

IF(NPUFF(MUNGRPI  .LE  .0)00  TO  250 

210  CONTINUE 

NPtl-NPUFF(MUNGRP) 

FLOCX(MUN)-XCENT(MUNGRP.  I  )+XB(MUN) 
FL OC Y ( MUN ) - YB ( MUN ) 

• FLOCZ(MUN)-ZB(MUN) 
FSI2X(MUN)-SIOX(MUNGRP, I ) 
FSIZY(MUN)-SIOY(MUNGRP.  1  ) 
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PSIZZ(MUN)-SlOZ(MUNORP,  I ) 

^i.OCX(MUN)  •XCCNTfHUNORP.NPM)  ♦XB(MUN» 

»^OCY(MUN)*YB(MUNI 

rCOCZrMUN)-ZB(MUNI 

NS  I ZX  f  MUN ) -S I  OX « MUN  GRP ,  NPM I 

NS  I ZY « MUN I -S I OY  ( MUNORP ,  NPM ) 

NS  I ZZ  ( MUN )  -  S I OZ  ( MUN  ORP ,  NPN ) 

CO  TO  P9g 

P90  CONTINUE 

NPUFFfMUNORP»-0 
FLaX(MUN)-0.0 
FLaY«MUN)-0.0 
FLaZlMUN»-0.0 
rsizx(MUN)"0.o 
FSIZY«MUN)-0.0 
rSIZZ(NUN)-0.0 
OCX  I  NUN) -0.0 
M.0CY(NUN)-0.0 
M.0CZrMlN)-0.0 
NSIZX(NUN)-0.0 
NSIZY(NUN)^0.0 
NSIZZ(MUN)-0.0 

P9g  CONTINUE 
RETURN 
END 
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SUBROUTIIC  DCNSTY(XP,YP,ZP.CON) 

COrfW  XBf  100) ,  YBf  I  00)  ,ZB(  100)  ,TB(  100) , 

COMMON  /AOEINO/  tAGEdOO) 

COMMON  nCSCR/  XCENT (30 . 150) 

COMMON  ’/HCSCR/  S30X  (30, 150)  .SIOY(30, 150  )  ,SIOZ(30, 150) 
COMMON  /HCSCR/  TACT T(30, 150) ,NPirr(30) 

CON«0 .0 

no  3QQ  MUNil .NBURST 


MUNORP-lAOEfMUN) 
NPM«NPUFF(MUN0RPI 
lF(NPM.LE.0)OO  TO  399 
X«XP-XB(MUN) 

Y-YP-YB(MUN) 

Z«ZP-ZB(MUN) 

CALL  OENM(MUNGRP,X.  Y.Z.C0t#1) 
C0N-C0N^C0»#1 
399  CONTINUE 

RETURN 

END 
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5UBR0UTI^€:  OCNMtMUNP,XP,YP,ZP,CON» 


COMMON  /HC5CR/  XCENTcao, |50) 

COMMON  /HC5CR/  SIGX J20. 150) ,SlOYiao, 150 > ,SIGZ(20, 150) 
COMMON  /VCSCR/  FACTTC?0,I50),NPUrFC80) 


COMMON  /POINT/  X.Y.Z 


X«XP 

Y«YP 

Z*ZP 

MUN-MUNP 


CON«0.0 

NPM-NPUTFCMtW) 

IFCNPM.LE.OIGO  TO  399 

IF(X.LE.XCENT(MUN.NPM)“4.0'SIGX(MUN,NPM)  )GO  TO  399 
IF  IX.GE.XCENTIMUN. I ) ♦H.O'SIGXIMUN, 1 ) )GO  10  399 


ir(NPM.GE.5)G0  TO  350 


DO  310  I“l .NPM 
CALL  C3ENPF(MUN.  I.TERM) 
CON-CON+TERM 
310  CONTINUE 
GO  TO  399 


350  CONTINUE 

DSTBPF-(XCENT(MUN. I  ) -XCENTIMUN.NPM) ) /FLOAT  (NPM- I  1 

NMID-(XCENT(MUN.  11  -  X) /DSTBPF+a. 0 

M-EFT-NMID 

NRIGHT-NMID-! 

IF  (^t-EFT.LE.O)^t-EFT-! 


360  CONTINUE 

IF(^t-EFT.GT.NPM)GO  TO  370 
CALL  DENPF(MUN,^t-EFT.TERM) 
CON-CON+TERM 
NLEFT-M-EFT+I 
IFCTERM.OE. 0.01 100  TO  360 

370  CONTINUE 

IF  (NRIOHT.OT.NPMINRIGHT-NPM 
375  CONTINUE 

IF(NRIOHT.LE.O)0O  TO  399 
CALL  DENPFIMUN.NRIGHT.TERM) 
CON-CON+TERM 
NRICHT-NRIGHT-I 
IFtTERM.GE.O.OllGO  TO  375 

399  CONTINUE 
RETURN 
END 
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5UBR0UTIIC  OENPFIMUf#».IP,CON) 

COMMON  /HCSCR/  XCENT<80, 150) 

COMMON  /HCSCR/  SIGX  <80. 150)  ,SIGYI80. 150  1  ,5102180. 1501 
COMMON  /HCSCR/  FACT  T<80. 150)  ,t#*Urri80) 

COMMON  /POINT/  X.Y.Z 

I-IP 

MUN-MUNP 

CON-0.0 

C0NN-<-.5*(  <<X-XCENT(MUN,I))/SIGX<MUN,  I  ))«»8 

•  ♦(Y/SI0Y(MUN,I))*«8 

•  ♦I2/SI02<MUN.I))«*8)) 

iriCONN.LE. -85.)  00  TO  350 
CON-FACTTIMUN.  I )  •EXP<CONN) 

350  CONTINUE 
RETURN 

END 
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SUBROUTINE  LCON  (POXP,POYP.POZP,l«P."*'P.UZP,TOrLNCI 
COMMON  X0( 1001 ,Y0f I  001 ,Z0C 1001 ,10(1001, N0URST 
COMMON  /A(3FIND/  lAOEdOOl 
COMMON  /HCSCR/  XCENT (20, 150.1 

COMMON  /HCSCR/  SIGX (20. 1501 ,SI0YI20. 150 1 ,SIGZ(20, 1501 
COMMON  /HCSCR/  FACT T(20, 1 501 .NPUTF (201 

UX-UXP 

UY-UYP 

UZ-UZP 

TOTLNC-0.0 

DO  ‘♦99  HJN-l  ,NBURST 

MUNGRP-IA0C(MUN1 
NPM-NPUrF(MUNGRPl 
IF(NPM.LE:.0)G0  to  H99 

POX«POXP-X0(MUN) 

POY-POYP'Y0(MUN) 

POZ-POZP“Z0(MUN1 

CALL  LCONM(MUNORP,POX,POY.POZ,UX,UY,UZ,  TERM) 
TOTLNC-TOTLNC+TERM 

*499  CONTINUE 

RETURN 

END 
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SUeROUTINE  LCONM  IMUNP,POXP,POYP,PO/P.DXP,DYt’,DZP,TOTLNC) 
COMMON  /HCSCR/  XCENUPO.ISO) 

COMMON  /HCSCR/  SI GX (PO , 1 50) .SIGYCPO . 1 50 ) .S IGZCPO , 1 50 J 
COMMON  /HCSCR/  FACT T ( PO . 1 50 ) , NRJFr C ao ) 

COMMON  /Lir«:/  POX,POY,POZ.DX.DY,DZ 

POX'POXP 

POY^POYP 

POZ-POZP 

DX-DXP 

DY*DYP 

DZ-DZP 

MUN-MUNP 

NPM*NPUFF(MUN) 

TOTCL-0.0 

ir(NPM.LE.O)GO  TO  H90 
VL*XCENT(MUN,NPM) 

VU-XCENTtMUN.I) 

C  PUFF  CHOOSING  LOGIC 

C  DETERMirC  INTERSECTION  WITH  Y.Z  PLArC 

IF  CDY.EQ. 0.0)00  TO  HIO 
S--POY/DY 
XINTER«POX+S*DX 
IF(XINTER.LE.VL)GO  TO  HIO 
IFIXINTER.OE. VU)00  TO  HIO 
DSTBPF-  ( VU-VL )  /FLDA  T  ( NPM- 1 ) 

NM I D-  ( VU-X I NTER )  /DS  TBPF+2 . 0 
KLEFT-NMID 
NRIGHT-NMID-I 
IF  (M_EFT.LE.0)KLEFT-1 
GO  TO  H70 
HIO  CONTINUE 
NMID«NPM/8 

CALL  LCONPF(MUN,NMID,TERM) 

IF(TERM.LE.O<OI )00  TO  H20 
KLEFT-NMID+I 
NRIGHT-NMID 
GO  TO  H70 
waO  CONTINUE 
KLEFT-I 

CALL  LCONPTIMUN.KLEFT.TERM) 

IFITERM.LE.O.OnOC  TO  H30 
NRIGHT-0 
GO  TO  H70 
H30  CONTINUE 
NRIGHT-NPM 

CALL  LCONPF(MUN.NRIGHT.TERM) 

IF  (TERM.LE.O.OI)  GOTO  H^O 
M_EFT-NPM+I 
GO  TO  H70 
HUO  CONTINUE 
TOTCL-0.0 
GOTO  H90 

H70  CONTINUE 

ir(KLEFT.GT.NPM)00  TO  HBO 

C  COMPUTE  LI^C  CONCENTRATION  CONTRIBUTION 

I-M_EFT 

CALL  LCONPF(MUN. I .TERM) 

TOTCL-TOTCL^TERM 

kleft-nleft^i 

IF(TERM.GE.O.OI )G0  TO  H70 
HBO  CONTINUE 

IF  (NRIGHT.GT.NPM)NRIGHT*NPM 
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‘♦95  CONTINUE 

IFfNRIOHT.LE.OlOO  TO  ‘♦90 
l-NRIOHT 

CALL  LCOf^JMUN,  I  .TERM) 
TOTCL-TOTCL+TERM 
NR10HT-NRI0HT-1 
iriTERM.GE.O.OI )GO  TO  405 
490  CONTINUE 

TOTLNC-TOTCL 

RETURN 

END 
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5U0ROUTI^C  LCONPT  (  MUNP, IP.PTLNCI 
COMMON  /HCSCR/  XCENKSO,  1501 

COfW^  /HCSCR/  5I0X  120, 1501  ,5IOY(20.  J50  J  ,5!0Z(20, 1501 
COf^lON  /HCSCR/  FACTT(20.I50)  ,NPUrr(201 

COMMON  /Lt^C/  POX,POY.POZ,OX,DY,DZ 

I-IP 

MUN-MUNP 

PFLNC-O.O 

A-(DX/5I0X(MUN, I  1) • *2^ 1 DY/5I OYC MUN, I  1 J • •2^CDZ/5I6Z(MUN, I ) 1 • *2 
B»DX»1P0X-XCENT(MUN,I) l/5tGX(MUN, I»**2* 

•  DY-(POY  )/5IGY(MUN,I)»*2* 

•  DZMPOZ  l/SidZCMUN,  1 1--2 

5--B/A 

PMX-P0X*5*DX 

PMY-P0Y*5»DY 

PMZ-P0Z*5»DZ 

C-( (PMX-XCENT(MUN, I  n/5I0X(MUN. I ) 1 ••2^ 

•  ((PMY  J/5I0Y(MUN.|IM«2+ 

•  ((PMZ  )/5I0Z(MUN,I»)**2 

CC-(-.5»C) 

ir(CC.LE.-25. )  GO  TO  399 

PFLNC-FACTT  (MUN,  I  >  •  2.5066283»EXP(CC  1  /SORT  (A) 

399  CONTINUE 
RETURN 

END 
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